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THEORY OF TRANSPORT IN LINEAR BIOLOGICAL SYSTEMS: 
I. FUNDAMENTAL INTEGRAL EQUATION 


JOHN L. STEPHENSON 
LABORATORY OF TECHNICAL DEVELOPMENT 


NATIONAL HEART INSTITUTE 
BETHESDA, MARYLAND 


The conditions under which the output, ¥p, (4), of a biological system is 
related to the input, y,(¢), by an integral equation of the type yp) = 


t 
[ Yq (@) w(t — @)d@, where w(t) is a transport function characteristic of 
0 


the system, are analyzed in detail. Methods of solving this type of inte- 
gral equation are briefly discussed. The theory is then applied to prob- 
lems in tracer kinetics in which input and output are sums of exponentials, 
and explicit formulae, which are applicable whether or not the pool is uni- 
formly mixed, are derived for ‘‘turnover time’’ and ‘‘pool’’ size. 


1. Introduction. Biologically the primary problem in which we 
are interested is the unique characterization of the transport of 
molecules from one part of a biological system to another. A mole- 
cule of a particular species may leave a particular volume element 
either by diffusion or convection or by some chemical reaction. 
The fundamental idea underlying the analysis of this problem is 
that of continuity; i. e., the rate of accumulation of particles of a 
particular species in a particular volume element equals the rate 
at which they are entering minus the rate at which they are leav- 
ing. As a corollary we have that the total number of particles in 
a particular volume element is equal to the summation with respect 
to time of the net flux of particles. Ordinarily it is customary to 
separate the flux term into a spatial term and a source term, de- 
pending on chemical reactions, within the particular volume ele- 
ment being considered. When this is done we obtain the classical 
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differential equation 
— div J, (r,t) + 8, (r,t) = 0¢,(r, 0)/dt, (1) 


where J, (r, t) is the total vector flux of particles of the #th species 
at position r and time ¢, s,(r,¢) is the ne¢ production per unit vol- 
ume from chemical reactions, and c, (r, ¢) is the concentration. 

Solutions of the above equation for biological problems have 
been considered by many workers (Rashevsky, 1948), and the litera- 
ture on mathematically related problems in diffusion and heat flow 
is tremendous. However, a useful closed solution of the equation 
depends on a reasonable geometry and a reasonable source func- 
tion. Frequently neither of these conditions is satisfied in a bio- 
logical system. Fortunately, for many purposes, solutions which 
are continuous functions of both the time and space variables are 
not necessary. For example, the problem might be to analyze the 
release of a metabolite from some particular organ of the mammalian 
body and its transport to another. Here, we are interested in total 
quantities of the metabolite in the two organs and the fluxes into 
and out of them; we are not at all interested (at least as a primary 
problem) in the distribution of the metabolite within the organs or 
within the blood stream. Problems of this kind have led to the 
abandonment of solutions in terms of the spatial variables and to 
various approximation methods. 

By far the most widely used approximation has been to divide 
the biological system into compartments. For purposes of analysis 
these compartments are considered to be uniformly mixed with con- 
stant fractional turnover rates. This approximation has consider- 
able anatomical and physiological justification for some problems. 
For example, we clearly have the intracellular and extracellular 
fluid compartments and if, for practical purposes, the only barrier 
to exchange of some substance between them is the cell membrane, 
a two-compartment approximation of this type should work very 
well. In practice it frequently does, and experimental data can be 
analyzed in terms of the expected two or three exponential curves. 
In other problems, while anatomically we still have the two or three 
compartments, the analysis of transport between them in terms of 
fractional turnover between uniformly mixed compartments is not 
even remotely correct. 

The mathematical nature of the difficulty becomes apparent when 
we derive the approximation of n uniformly mixed compartments 
from the fundamental equation (1). As a first step in this deriva- 
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tion we integrate both sides of equation (1) over a finite volume to 


give 
-{ div J, dV, -{ s, dV, -/ dc, /dtdV, . (2) 
V; V; V 


i t i 


If we utilize the divergence theorem, we can replace the volume in- 
tegral of the divergence of J, by the surface integral of J, over the 
surface S, bounding the volume Vs Doing this and inverting the 
order of differentiation and integration on the right-hand side of 
equation (2), we obtain 


= J, -da+ | s,dV,=dq,,/dt, (3) 
ye ae 


u 


where q,, is the total amount of material of the &th species in the 


zth compartment. (It will be recalled that -{ J,:do is the net 
Sh 
inward flux.) 

As it stands, equation (3) is no more amenable to solution than 
is equation (1). If, however, the integrals on the left-hand side 
can be replaced by more manageable expressions, an exact solu- 
tion may be possible. The simplest approximation is to replace 
them by linear terms in the quantities q,, the so-called uniformly 
stirred multi-compartment approximation. Here, we assume that 
the flux of particles of the &th species from the jth into the zth 
compartment is given by a product ki, Yk? where kj, is called 


the fractional exchange coefficient or the fractional ane rate. 


We then have 
-{ Jdo=)) ki Tjk (4) 
S J 


i 


where k,,, = =Do, . Approximation of the source integral 
Md meee si 
by a similar sum gives 


He 5a, =D) Fear (5) 
t 


where k. ikl is the fractional conversion per unit time of the /th to 
the kth species in the ith compartment, and where k;,, = 


“ae k, lk Thus we obtain the linear first-order differential 
Ldl¢k) ’ 
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equation for q,, 
dq,,/dt =), Kaye Vin + 2 Peete , (6) 
J 


and similar equations, nm in number, over the n compartments and 
the m species. Such linear systems of equations have been dis- 
cussed by numerous authors in connection with the transport of 
radioactive label (Robertson, 1957). They are also the subject of 
a large classical physical and mathematical literature. In a wide 
class of problems, the application of these equations has been 
very successful. However, in certain problems the above approxi- 
mation of the surface flux integral is not even approximately valid. 
One such is that of exchange between two compartments connected 
by a circulatory system, e. g., the right heart and the left heart. 
To handle this type of problem, a transport or transfer function, 
w(t — @), can be introduced which when multiplied by AZ¢ defines 
the probability that a particle which enters the first compartment at 
time w enters the second, for the first time, during the time in- 
terval Az. For the simple case of a circulatory bed this leads 
(Stephenson, 1948, 1958; Sheppard, 1954, 1957) to an integral 
equation, namely, 


r(0 = [ y.(w)w(- 0) de, (7) 


0 


where y(t) is the flux of particles entering the first compartment 
and y,(¢) is the flux of particles entering the second compartment. 
Integral equations of the type of (7), frequently known as 
Volterra’s equation (Whittaker and Watson, 1943, p. 221) arise in a 
variety of problems: for example, ‘‘aging’’ processes, (Shemin and 
Rittenberg, 1946; Ottensen, 1955; Bateman, cit., Titchmarsh, 1937) 
and in “‘linear filter’’ theory (see James, Nichols and Philips, 
1947, chap. 2). A modification of this equation has been applied 
to metabolizing systems by H. Branson (1946, 1947). Branson’s 
particular formulation has been criticized by others (Wijsman, 1953; 
Hearon, 1953). As we will show in detail below, (7) is the charac- 
teristic equation of any system whose output summates linearly 
both with respect to time and to input. It is the purpose of this 
and following papers (Stephenson, in preparation) to develop a 
systematic mathematical procedure for applying this probabilistic 
or integral equation method of analysis to transport in biological 
systems. In this paper we will first formulate the prototype of 
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equation (7) for a two-flux system, analyzing in detail the neces- 
sary conditions of linearity. We will then discuss general methods 
of solution briefly and will obtain explicit closed solutions for 
some relatively simple (but frequently occurring) problems. In fol- 
lowing papers (idid.) we will extend the analysis to multi-flux sys- 
tems and will consider some general problems in the analysis of 
data and in the construction of models in systems with incomplete 
information. 


2. The fundamental equation. The prototype of the problem we 
are considering is shown in Figure 1, in which we have shown an 


w 

ats 

FIGURE 1. Influx, vs. (¢), and efflux, Y, (¢), to a region of arbitrary 
geometry. 


irregularly shaped region in space which is bounded by an impervi- 
ous surface, except for a surface of entry, S,, and a surface of 
efflux, S,. We will suppose that a large number of molecules M, 
pass at time ¢, through S,. Initially, to simplify analysis, we will 
suppose that eventually al/ of these molecules leave through S,. 
We will also suppose that it is possible to record the molecules 
which leave the system during the various time intervals At, = ¢, - 
t,_, after their entry, and that once the molecules leave the sys- 
tem they do not re-enter it. The result of such an experiment could 
be represented by the series of rectangles shown in Figure 2, the 
number of particles leaving the system between ¢,_, and ¢, being 


N,=M,w,At;, (8) 
this equation defining w,. By such an experiment we could define 


a ‘‘function’? w(M,,¢,,¢—¢,), Shown by the smooth curve in Fig- 
ure 2, such that the number of particles which leave the system 
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to t, te ts ty ts th-2 thet ie 
{+ —_———q@q“~ 


FIGURE 2. Transport function w(t). The rectangles are proportional in 
area to number of particles in efflux during indicated time intervals after 
the sudden introduction of material as a bolus at time ¢). 


during the time interval Az would be given by 


N, = My (Mos tort — to) At. (9) 


Such a ‘‘function,’’ although it can be defined for a single ex- 


periment provided that the operation leading to the computation of 
the numbers w, can be carried out, is useless unless it is in some 
way characteristic (i. e., mathematically invariant) for the system. 
As we have indicated in our notation, the numbers w, we obtain 
as the result of an experiment may depend not only on the time 
difference ¢—¢, between entry and exit, but possibly on the num- 
ber of particles which enter the system and on the time of entry. 
Other parameters which we have not indicated, such as the dis- 
tribution of the particles in space and momentum, may also affect 
the numbers w,. At worst, we may be able to discern no relation 
between influx and efflux other than conservation of matter. A 
priort, for an arbitrary system we have no way of deciding upon 
what variables the response of the system is going to depend. 
However, we can approach the problem experimentally by varying 
the input with respect to time and amount. If in repeated experi- 
ments (with sufficient time between experiments for all the mole- 
cules to leave the system) the numbers w, are invariant with re- 
spect to the exact time and amount of the input, they will define a 
‘‘normal’’ response or transport function such that 


N,=M,w(t-t,) At. (10) 
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This response function is in a sense characteristic for the system; 
however, we must consider the additional complication of an input 
which is not a discrete group of particles at a single instant of 
time but which is distributed in time. Suppose, for example, the 
input consists of M, particles at time ¢, and M, particles at later 
time ¢,. In general the response of the system to the second group 
of particles may be modified by the prior introduction of the first 
group, but in certain systems, the system will have the same normal 
response to the second influx as to the first; that is, of the parti- 
cles in the efflux during the time interval At, those contributed by 
the first input will be M, w(t — t,)At, and those contributed by the 
second input will be M, w(t —¢,) At, where w is the same function 
for both inputs. For the total efflux due to both inputs we have 


N,=M,w(t—t,)At+M,w(t-t,) At. (11) 


As an obvious generalization of equation (11), for certain systems 
we will find that for an input distributed in time, with amount M, in- 
troduced at time ty the efflux will be given by 


N, is od RSE (12) 


where the efflux for each input is the same function of ¢ -¢,. A 
system whose response satisfies (12) is said to be linear. 

From equation (12) we easily derive our integral equation. If 
we first define an outgoing flux function y, (¢) by 


N, =y,(é)At, (13) 


we obtain 


y(t) At =)" Mpwt-t) At. (14) 
J 


Since the A¢ on each side of equation (14) refers to the same time 
interval, they both may be cancelled to give 


y, (2) =)! M,wt-t,)- (15) 


Equation (15) holds for an arbitrarily large number of discrete in- 
puts. If we define an input function y, (¢,) by 


My =ye(t)At;, | (16) 
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where M. is now the total number of particles which enter the sys- 
tem during the time interval AZ,, and substitute (16) into (15), we 
obtain 


yp) =)" valt)we-t)At,. (17) 
J 


By dividing the interval between ¢, and ¢ into arbitrarily small in- 
tervals (we assume there is no entrant flux before time ty), we 
have as a limit, setting ¢, = 0, 


t 
ie [ y, (t,) w(t - t,)dt,. (18) 
0 
The simple transformation = ¢ — t gives the alternative equation 


n= [oy (¢ a) w(o) da. (19) 
0 


Before discussing the solution and application of this integral 
equation, we will consider further the system of linear algebraic 
equations from which it has been derived as a limiting form. This 
will provide insight into the difficulties which arise when the sys- 
tem is nonlinear. If as we have assumed above, the interval from 
t, to ¢ is divided into n sub-intervals and if N,. is the number of 
particles in the outgoing flux during the time interval Az, which 
are contributed by particles from the flux M, which enter the sys- 
tem during the time interval At;, j Si, we can clearly define a 
number W sj by the equation 


N,,=M.w,,. (20) 


ty oS 

Furthermore, if by some kind of distinctive labeling we could in 
some way identify the molecules which entered the system during a 
given time interval A¢., we could actually measure the numbers 
Ww; 5 for a particular set of M.. However, as we pointed out above, 
for an arbitrary system every different set of M. may give a dif- 
ferent set of © ,;. OF even the same set of M. used in two different 
experiments on ‘the same system may yield different numbers W;; 
Only for special systems will the w,, turn out to be invariant with 
respect to the input. 

In practice, we usually cannot identify the input by character- 
istic labeling for every time interval Az j and can only measure the 
total efflux and influx of material. The total efflux during the time 
interval At;, N;, is obtained by summing over all previous influx, 
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N, = N= Dd) Ms (CS Gee) (21) 


If we take the sets of N, and M, obtained for a particular experi- 
ment and assume they are related by such a set of linear equa- 
tions, unless we place further restrictions on the w,,, we have n 
equations in (n + 1)n/2 unknowns. Naturally this system does not 
have a unique solution. (Under certain circumstances, for example, 
if the w,,; were independent of the M, and depended only on the 
time of input, we could get a unique solution by performing a num- 
ber of experiments greater than (nm + 1)/2, each with a linearly in- 
dependent set of inputs M ;-) However, if we place the further re- 
striction on the ws that the system is linear as defined above, 


that is, Wi, = W,, ifi-j =r—s, so that 
ah AE pee RA Dy 7 = 9; 
OSG cP eR Larne Sh oe oe he nale ae o (22) 
own = Miami = Ma-1» 
bad oe = Wns 


we obtain a system of n equations in the n unknowns: w,,w,...W 
specifically: 


n? 


2 a> 1 q = (23) 
N. =M w,+Mi- Waters +M,w,. 


If the determinant of the M, is non-zero (i.e., M, #0), the system 
(23) can be solved uniquely.* 


*The author is indebted to J. Z. Hearon for pointing out that the 
equivalent linear set for 


t t 
Vp -{ y_(¢= 0) (ada = [ y(w) w(t-w)do 
ty) 0 


is triangular because the integration is from zero to ¢. But given tri- 
angularity the restrictions (22) result because the arguments t-@ and 
@ can be interchanged in a convolution integral. Otherwise, e. g., if 
w is a function of @ and t—@, this is not permitted. Thus convolution 
(whose basis is linearity or superposition) insures triangularity and 
equality on the diagonals of the determinant | M |. 
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The pitfall in this solution, which is also the pitfall in any solu- 
tion of the corresponding integral equations (18) and (19), cannot 
be over-emphasized. This is that any single experiment which 
leads to a set of N, and M, will give a solution for the w,, if the 
assumptions embodied in (22) are made. However, unless the as- 
sumption is valid, the w, are in no way characteristic of the sys- 
tem. In practice one must either have some prior knowledge of the 
system which makes the assumption of linearity reasonable or 
else compute the w, from a number of linearly independent sets of 
M.. Usually some prior knowledge of the system is assumed: for 
example, the transport of trace amounts of radioactive label in 
steady-state biological systems is assumed to be a linear process; 
or to draw examples from other areas, the response of an electrical 
filter composed of linear circuit elements or the passage of in- 
dividuals through a population pool is taken to be linear. If the 
system is truly unknown, an exhaustive approach to determine its 
characteristic response under all distributions of input would be 
impractical. For example, if the period from ¢, to ¢ is divided into 
ten sub-intervals and if the M, for each interval can assume inte- 
gral values 1 through 10, there will be 10!° possible sets of M.. 
For the same reason the chance that two different sets of M, will 
give the same values for the w,, unless the w, are characteristic 
for the system, is extremely small, and one would ordinarily be 
satisfied that the w, were characteristic of the system if a few in- 
dependent sets of M. gave nearly the same values for them. 

The derivation and the limitations of the above approach to 
transport problems have been gone into in such detail both be- 
cause of the controversy which has arisen over its application to 
metabolizing systems (Branson, 1946, 1947; Wijsman, 1953; Hearon, 
1953), where in general the conditions of linearity are not satisfied, 
and because nowhere have the concepts underlying its application 
to biological systems been fully developed. 


3. Solution of the integral equation. In this section we will 
briefly outline the standard methods of solving (Mikhlin, 1957; 
Whittaker and Watson, 1943, chap. 11) the integral equation 


a (= [y= o)wlo)do. (24) 
0 


As is obvious from our derivation, one method is to approximate it 
with a finite system of linear algebraic equations. Such an ap- 
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proximation is practically important in numerical solutions with 
computers and also underlies some of the analytical methods. 

By differentiation with respect to ¢, equation (24) can be trans- 
formed into an equation of Fredholm’s type (idid., p. 221): 


yi () <y7, (Otol) + fe SEN (25) 


0) 


A solution of (25) is a solution of (24) provided y, (0) =0. One 
method of solving (25)—not of great practical importance but of 
considerable interest—is that of successive substitutions (idid., 


p. 221). Setting y, “(¢)/y, (0) = f(¢), and Yq (- )/y,(0) = 9(¢-@), 
we have 


w(t) = f() - [ 9@- 0) w(w) de. (26) 


0 


Successive substitutions give: 


v0=1- [ 9e-o)fe)de, 
0 


w, (2) -1)- [ e(t- 0), (ado, 
0 


(27) 


w(t) = (2) - [ee -—o)w,_,(o)do. 
0 


This sequence of functions can be shown to converge to a limit for 
finite values of ¢, and by direct substitution it is seen to be a so- 
lution of (26). 

Ordinarily, however, an analytic solution of equation (24), if it 
exists, is most simply obtained by the use of Laplace transforms 
(Churchill, 1958; Sneddon, 1951; Titchmarsh, 1937). The use of 
transforms also greatly simplifies the discussion of multi-compart- 
ment systems. By definition, the Laplace transform of a function 


f(t) is: 
Lf() = F(p) = i ” #(t) Pde, (28) 
10) 


where p is a number whose real part is positive and large enough 
to make the integral convergent. (In these eauations and else- 
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where in the paper, functions are denoted by lower case letters 
and the transforms of functions by capital letters, unless it would 
lead to confusion; then the transform is indicated by a preceding 
L.) The fundamental properties of transforms of which we will 
make frequent use are 


Lf’ (t) = pF - f(9), (29) 
L [roar Eire (30) 
0 
and for two functions /, (¢) and f, (¢), equal to 0 for ¢ < 0, 
Lf ne ~ 0) f, (0) do. = F Fy; (31) 
0 
also, formally,* 
pos -[ f(t) dt, (32) 
0 
-(dF/dpl,- , -[" tf (t) dt. (33) 
0 


Application of (31) to the integral equation (24) gives: 


W =I,/T,. (34) 


From equation (26) we obtain, by taking transforms of both sides 
and solving for the transform of w, 


W=F/(1+ 9). (35) 
Provided that y,(0) = 0, this reduces to (34). Aliso, if ®(p) is less 


than unity in absolute value, (35) can be expanded to give 
Wm F(1- 6407 =.) G1 I (36) 


Equation (36) can also be obtained directly from the iterated solu- 
tion (27). 


4. Applications. As a first example of the application of the 
formulae derived in the preceding section we will consider a sys- 


*Although for actual convergence of integrals, restrictions must be 
placed on the class of admissible functions. 
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tem in which the influx and efflux are related by an equation of 
the type of (24) and are both sums of exponentials: 


valde exds ent, (37) 
—bt 

y,(¢) Sey j (38) 
The transforms of (37) and (38) are 

i; =) A/(p +4,), (39) 

3 =2, 8/0 +5). (40) 
From (39) and (40) we obtain 

Pr, =P,/IL(@ +4), (41) 

DT, = P,/JT,( + 4;), (42) 


where if there are n terms in y,(¢), P, is a polynomial in p of de- 
green —l,i.e., 


Beg Pr et nd os PP ase HP 


a | amd 


(43) 


0 pb] 


where 


Pa. =) Ae (44) 


= ’ a 45) 
me 7, A | ( 
PA =)" ‘t4, = (all & - 1 products of different a, q&i)]. (46) 

t 


Similarly, if y,(¢) has m terms, P, is a polynomial in p of degree 

m-—1 or m—2 (usually m — 2 since ordinarily > iB = 0 because 
j 

y,(0) = 0). Substituting (41) and (42) into (34), we obtain 
4 P, JI, + a) 


S ’ (47) 
= HT (p + b,) 
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or denoting the n - 1 roots of P, by 7,, 
P, IT,” +4) 


feel eile ee eee Se (48) 
IT,(? + 1;) I1,(@ + 4) 


If the fr and b, are all distinct and ) | B, =0, equation (48) can be 
j 
reduced to the partial fraction form 


W = pee C,/(p + ¢,); (49) 


where the c, are the n + m—1 quantities r, and b;. Equation (49) 
can be inverted to give 


w(t) aes ie De (50) 


Frequently we can get useful results without going through a 
complete solution of the problem. For example, we have from (34) 
and (32) 


ik w(t)dt = ees “[y7 8/8) | [yo 4/1] ‘ (51) 


In the simple conservative system that we have been considering 


we should have | w(t)dt =1; however, in multi-flux systems this 
0 


is not necessarily true. 

Another useful formula is that for computing the average time a 
particle Jeaving the system has spent in the system. We have the 
general formula 


2 S dW 
t t)dt = rs ¢ eae fs (eae 
ie w (t) I t w(t) dt el (52) 


where i w(t)dt £1. From (34) we obtain 
0 


dW [dU ,/dp a /d 
Wee eis |, (53) 


dp ie r 


a 
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whence 
i ty, (t)dt | ty (t)dt 
0 
2 eee career (54) 
if y,(t) dt | y,(t) dt 
0 0 
or, 
(ied pat (55) 


where t, is the weighted average time a particle leaves the system 
and ¢, is the weighted average time of entry. Equation (55) is ob- 
viously valid, but its derivation from (34) illustrates the manner in 
which relations in Laplace transforms go over into integrated fluxes 
and probabilities for p = 0. It will be observed that (54) is valid 
regardless of the functional form of y,(¢) and y,(¢). In the particu- 
lar case we have been considering in which the fluxes are sums of 
exponentials, either by differentiation of the transforms or direct 
application of (54), we obtain 


Mei EeS Ada? 


oS, Se eee ee ee (56) 
Se SRY es SOT Aya: 


When the efflux is partitioned into several fractions with corre- 
sponding partition of the transfer function, i.e., 


y,() = Zy,(), (57) 
w(t) = % w, (2) : (58) 
f v@aezf w,(t)dt = 2W,(0) = 1; (59) 
0 0 
we have 
bin > W, (0), -t, : (60) 


Also in a steady state conservative system which has a steady 
influx of material at constant rate m, starting at zero time, the 
total quantity of material in the system at time ¢ is: 


t Q 
(= me- | ao mw(o)da (61) 
0 ) 
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-m[e- [dof wtordo. (62) 
0 a) 


As t — «, equation (55) can be shown (Stephenson, 1958) to re- 


duce to 


5) 


q(~) = mt,,- (63) 


Thus we are led to a straightforward method of computing total 
‘pool size’? and ‘‘turnover time’’ directly from experimental data 
without reference to any particular compartmental model. It should 
be noted that ¢| from (54) is independent of the time course of in- 
troduction of labeled material into the system; m of course corre- 
sponds to the total flux of unlabeled material through the system. 

Although some of the above results can be obtained without the 
use of transforms and transfer functions, this section has illus- 
trated their application to a familiar biological problem. It should 
be emphasized that the primary object is to obtain the transfer 
function w(t) from the fluxes and not the opposite. Practically 
such a function derived under certain experimental conditions can 
be used to compute fluxes not directly measurable under other 
circumstances. One problem in which this approach has been found 
very useful is in comparing the metabolism of C** labeled un- 
esterified fatty acid and C’‘ labeled tri-glyceride fatty acid (Stephen- 
son and Fredrickson, unpublished). 


The author thanks Mr. Arnold Jones for his general assistance 
and Dr. John Z. Hearon for his very helpful comments on this 


paper. 
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A review of the literature of aortic elasticity data leads to the con- 
clusion that the pressure-volume curves of the human aorta can be reas- 
onably well approximated by parabolas, while a straight line approxima- 
tion proves inadequate. By using the above relationship and assuming 
flow through the peripheral vascular system to be directly proportional to 
the pressure difference across it, two equations are derived which relate 
systolic and diastolic blood pressures, cardiac output, pulse rate, pe- 
ripheral resistance, and two constants of aortic elasticity. Pulsewave 
transmission phenomena are ignored. These two equations can be solved 
simultaneously for any two variables in terms of the others. When values 
computed from these equations are compared with experiment, the errors 
are uSually less than 5 per cent, provided errors due to inaccurate meas- 
urements of the aortic elasticity constants are avoided. 


Introduction. Progress toward a complete analysis of the mathe- 
matical physics of the human systemic arterial system has been 
made along two quite different pathways. The first approach con- 
sists of an analysis of pulsewave transmission in the aorta. This 
has been done in a rigorous fashion starting with the fundamental 
equations of hydrodynamics and using as boundary conditions the 
experimentally determined data on aortic elasticity (Wormersley, 
1957, 1958). Because of the mathematical difficulties inherent in 
this type of analysis, it has not yet been carried far enough to 
yield results of clinical usefulness. By a much less rigorous anal- 
ysis, O. Frank (1930) derived an equation for the volume of blood 
stored in the aorta with each heart beat in terms of aortic pulse- 
wave transmission data. K. Wetzler and A. Béger (1939, p. 474) 
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assume that this volume is equal to % of stroke volume and thus 
calculate the latter clinically useful quantity. R. Ronniger (1954~ 
55) handles the problem in a more rigorous way, using complex 
number theory. 

The other approach to the mathematical physics of the arterial 
system is a study of the system as a whole, assuming the aorta to 
be an elastic reservoir which expands and contracts simultaneously 
in all its parts. Equations describing arterial blood pressure as a 
function of time were derived by Frank (1899), using the assump- 
tion that dP/dV = constant (P = pressure, V = volume) for the aorta. 
Wetzler (1938) takes Frank’s (1899) equation for aortic pressure 
vs. time during systole (where flow from heart to aorta is assumed 
constant over time) and Frank’s (1899) pressure vs. time equation 
for diastole, and assumes an equilibrium condition to exist as fol- 
lows: At the beginning of systole, let aortic pressure =P ,. During 
systole, there is constant flow from heart to aorta causing aortic 
pressure to rise to a maximum (P_) at the end of systole when the 
aortic valve closes. Then the blood stored in the aorta gradually 
runs out through the smaller peripheral vessels until aortic pres- 
sure has declined to Lae at the end of diastole, when the aortic 
valve opens and the cycle is repeated. Under these conditions, 
Wetzler has derived essentially two simultaneous equations. 


Vink ts ae 
= ee ee 
a E jlT,E RR | (1) 
— . ah: 
Lae ed A Se 8 (2) 


P, = systolic pressure 

P , = diastolic pressure 

R = total peripheral resistance 

V_ =stroke volume = output of heart per beat 
T, = duration of systole (aortic valve open) 

T, = duration of diastole (aortic valve closed) 

E” =dP/dV (elasticity constant for aorta) 

Equations (1) and (2) can be solved simultaneously for either of 
any pair of variables as a function of all of the others. Wetzler 
(1938) eliminates P_ and solves for P, as a function of the other 
variables. He computes P, on various individual patients from 
measurements of the other variables and obtains reasonably good 
agreement of computed and experimentally determined values of 


a 


a 


a 
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P,. The aortic elasticity constant was obtained from pulsewave 
transmission data. Wetzler and Boger (1939, p. 530) give computed 
curves showing how P, varies when each of the other physical 
variables changes individually, with the other variables held con- 
stant. 

It was the opinion of the present author that equations such as 
the above, relating the physiologically important variables, might 
be useful in calculation of hemodynamic data in individual patients 
and might demonstrate clinically important relations between the 
physical variables of the human arterial system. Since the recent 
extensive data of J. W. Remington et al. (1948) has shown that the 
linear approximation of Frank (1899) describes poorly the pres- 
sure-volume curves of human aortas, it was decided to develop 
new equations using a parabolic approximation to the aortic pres- 
sure-volume curve. 

The over-all plan of this study is to determine from existing ex- 
perimental data the physical characteristics of each component of 
the human arterial system. Then, the system as a whole is ana- 
lyzed to determine the relation of each of the variables to all of 
the others. 


Physical Characteristics of the Components of the System 


A. The Aortic Pressure-Volume Curve. The experimental de- 
termination of the aortic pressure-volume (P-V) curves is most 
easily done if the aorta has been removed from the organism. 
Hence, many measurements have been made on the isolated aorta 
but very few on the aorta in the living animal. R. W. Lawton and 
L. C. Greene (1956) determined aortic P-V curves from measure- 
ments made on moving picture photographs of the movements of 
beads sewn on the outside of the abdominal aorta. These measure- 
ments showed the dog abdominal aorta to be considerably stiffer 
than did those of Remington et al. (1945) on the excised dog ab- 
dominal aorta. However, the method of determination of aortic di- 
ameters from photographs of beads is difficult and tedious, so that 
only a few sets of measurements were done. Since, in addition, 
the in vivo data was not compared with the P-V curves for the 
same aortic segments when excised, it would seem unwise to rely 
very heavily on these experiments for making any extrapolations 
to the intact living from the excised dead aorta. An attempt was 
made by L. N. Katz et al. (1947) to measure aortic P-V curves in 
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living dogs by using 0.1 second exposure X-ray pictures to de- 
termine aortic volume. Plots of pressure vs. volume thus obtained 
showed such wide scatter that the present author is reluctant to 
draw any conclusions from this work. After a review of the above 
literature, the author came to the conclusion that an adequate de- 
termination of the aortic P-V curve in the living animal had yet to 
be accomplished. 

The best available approximation to the above is the determina- 
tion of the P-V curves on an isolated dog aorta perfused with the 
blood of the same dog, as done by R. S. Alexander (1947). This 
experiment gives a P-V curve that is moderately stiffer than those 
of Katz et al. (1947) where the isolated aorta was perfused with 
isotonic saline solution, but the difference is not large enough to 
rule out the possibility that it might be due to a difference be- 
tween the individual dogs used. 

The experiments discussed in the preceding paragraphs suggest 
that the aorta in the intact animal is stiffer than when isolated 
from the animal, but these researches do not justify definite con- 
clusions. 

Another uncertainty in aortic P-V curve measurements was 
brought to light by the experiments of Lawton (1955) who showed 
that, for some conditions at least, the aortic elastic characteris- 
tics were quite different under static and dynamic conditions. 
However, it has been found that when dog aortas are subjected to 
physiological rates of stretch (Hamilton e¢ al., 1945), or when 
human aortas are expanded and contracted sinusoidally at a fre- 
quency and amplitude in the range of the fundamental frequency 
occurring under physiological conditions (Cope, 1958a), the P-V 
curves obtained are quite close to those obtained under static 
conditions. 

Hysteresis effects constitute another problem encountered in 
the measurement of aortic P-V curves. In the thoracic aorta of the 
living dog, R. F. Rushmer (1955) showed a P-V hysteresis loop in 
a direction indicating that the aorta is losing energy to the blood 
with each cycle of expansion and contraction, but hysteresis loops 
in the opposite direction were found by others (Lawton and Greene, 
1956; Alexander, 1947; Hamilton e¢ al., 1945). Because of these 
uncertainties, hysteresis effects will be ignored in this analysis. 

Of the various P-V measurements which have been made on iso- 
lated aortas, some used experimental set-ups allowing longitudinal 
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expansion while others did not. Longitudinal stretch of nine dif- 
ferent excised human descending thoracic aortas ranged from 4 to 
19 per cent when pressure was increased from 0 to 200 mm Hg 
(Cope, 1958b). In order to know which sets of measurements will 
best approximate the P-V curves of the aorta of the living organ- 
ism, it becomes necessary to know whether a change in pressure 
causes a change in the length of the aorta in situ. 

The human aorta in situ is apparently held under considerable 
longitudinal tension, since when it is cut loose from its attach- 
ments, it contracts longitudinally by 2 to 15 per cent (Remington, 
et al., 1948). It is therefore not surprising that in intact aortas in 
living dogs, changes in length usually produced volume changes of 
less than 25 per cent of those induced by diameter changes (Law- 
ton and Greene, 1956). It should also be mentioned that the length 
changes were 180 degrees out of phase with the pressure and di- 
ameter changes (Lawton and Greene, 1956). Based on this evi- 
dence, it would seem reasonable to suppose that aortic length 
changes are small in the intact aorta in sztu compared with the ex- 
cised aortic segment. 

There are available in the literature two reasonably compre- 
hensive studies of the P-V curves of excised human aortas as a 
function of the age of the individual (Remington ez al., 1948; 
Hallock and Benson, 1937), but the two used different methods of 
measurement and showed markedly different effects of age. P. 
Hallock and L. C. Benson (1937) inflated segments of descending 
thoracic aorta with isotonic NaCl solution and found large in- 
creases in both baseline volume and in stiffness with increasing 
patient age. Remington et al. (1948) measured force-length curves 
on a series of rings cut from different portions of the aorta and 
then computed a P-V curve for the whole aorta on the assumption 
that there was no longitudinal stretching. They found a progres- 
sive increase of baseline volume with increasing patient age but 
only slight changes in the shape of the P-V curve with age. 

Since, as mentioned before, longitudinal stretch probably does 
not occur to any great extent in the intact animal, the aortic P-V 
data of Remington e¢ al. (1948) will be used in the analysis. 

The data of Remington e¢ al. (1948) were published in the form 
of pressure vs. volume per meter? of body surface, so it was nec- 
essary to convert it to the form of pressure vs. volume for use in 
this analysis. Then all the volumes for each pressure were aver- 
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aged for each age group giving a pressure-volume curve for each 
age group as well as for each patient. The pressure-volume curve 
appeared to be roughly parabolic in shape with the axis of the 
parabola parallel to the V-axis as shown in Figure 1, so the equa- 
tion V = aP?+ bP +c (a, 6, c = constants) was fitted to the P-V 
data of each individual and each age group, substituting the pres- 
sures and volumes at 80, 140, and 170 mm Hg in the above equa- 
tion and solving simultaneously for a, 6, and c. Typical values 
fora, 6, and c are ~5.348 x 10° °,2.032, and -36.06, computed from 
the data of Remington e¢ al. (1948) for the mean values of the 20- 
30 year age group. The resulting curve fit is good from about 170 
to 50 mm Hg but is poor below 50 mm Hg. As an illustration of 
this for the means of the 18-33 year age group, at pressures of 20, 
50, 80, 110, 140, and 170 mm Hg, the measured aortic volumes are 
42, 65, 98, 125, 144, and 155 cc respectively. For the same pres- 
sures, the volumes computed from the approximating equation are 
2, 538, 93, 122, 144, and 154 cc. The maxima on most of the para- 
bolic approximations occur only a little above 170 mm Hg, after 


dv 
which aP * 0. Since above 170 mm Hg it has been found experi- 


mentally that the radius of the aorta remains almost constant 
(Lawton and Greene, 1956), the parabolic approximation of the P-V 
curve was used only for P< 170 mm Hg. At higher pressures the 


170 P 


FIGURE 1. The solid line shows the assumed shape of the aortic P-V 
curve. The parabolic approximation is used where P < 170 mm Hg. For 
P>170 mm Hg, a straight line approximation is used. The unused por- 
tion of the parabola is shown by the dotted line. 
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volume of the aorta was assumed to remain constant at the same 
volume as that given by the parabolic approximation at P = 170 mm 
Hg (Fig. 1). Because the parabolic approximation could not be 
used above 170 mm Hg, it was necessary to derive a separate set 
of equations to be used in this situation. 


B. Pressure-Flow Relations in the Peripheral Vessels. Classi- 
cal Windkessel theory, e.g., Frank (1899) and Yetzler (1938), as- 
sumes pressure to be directly proportional to flow in the peripheral 
vascular system. For the hindlimb of the dog, pressure-flow curves 
are linear (Pappenheimer and Maes, 1942; Whitaker and Winton, 
1933), but the pressure-axis intercept is positive and varies with 
the degree of vasoconstriction (Pappenheimer and Maes, 1942). 
The pressure-axis intercept, as well as the slope, also varies 
with the hematocrit of the perfusing blood (Whitaker and Winton, 
1933). With kidney, however, the pressure-flow relationship is 
quite different, with flow holding practically constant over a wide 
range of pressure (Kramer and Winton, 1939). For perfusion of the 
whole animal (spinal kittens), pressure-flow curves were nonlinear 
and changed shape over time as perfusion was continued (Folkow, 
1952), but shape changes as a function of vasoconstriction were 
not studied. Because adequate studies of pressure vs. flow for 
the whole animal as a function of vasoconstriction have not been 
made, this analysis will consider pressure to be directly propor- 
tional to flow as was done by Frank (1899). 

Most clinical investigators have assumed flow to be directly 
proportional to pressure difference with a constant of proportion- 
ality (R) called total peripheral resistance (an analogy to electri- 
cal resistance). The better determinations of Ff, e.g., Landowne 
et al. (1955), use time averages of arterial pressure divided by 
cardiac output using the direct Fick or the dye dilution methods. 
The commonly used calculation of mean arterial pressure as some 
linear function of systolic and diastolic pressures is subject to 
considerable error because the shape of the arterial pressure-time 
curve can change considerably as a function of the aortic P-V 
curve. It is also desirable to subtract an estimate of venous pres- 
sure from mean arterial pressure before dividing by cardiac output. 


C. Pulse Rate. Length of Systole as a Function of Pulse Rate. 
The standard method of pulse rate measurement is simple and di- 
rect and not greatly susceptible to error. However, the pulse rate 
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of any given individual varies greatly over time as a function of 
physical and mental activity (Addis, 1922; Boas and Goldschmidt, 
1932). 

With regard to the length of mechanical systole as a function of 
heart rate, A. A. Luisada (1953, p. 55) states that ‘‘mechanical 
systole is about 0.04-0.05 seconds shorter than electrical sys- 
tole.’? He also gives (p. 238) a table of values of the duration of 
electrical systole as measured by the Q-T interval of the electro- 
cardiogram as a function of heart rate. From these data were com- 
puted values of the durations of mechanical systole and diastole 
as a function of heart rate, as shown in Figure 2. The length of 
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FIGURE 2. Duration of systole and diastole as a function of heart 
rate. 


mechanical systole will be used here as the best available ap- 
proximation of the length of time that the aortic valve is open. 
Similar, but not identical, data is given by L. E. Praywara (1933). 
These figures should only be used with caution, since it has been 
shown in the heart-lung preparation that the duration of mechanical 
systole also depends on cardiac output and on the presence of 
sympathomimetic amines in the blood (Braunwald e¢ al., 1958). It 
is also reasonable to expect that abnormalities of the cardiac con- 
duction system might alter the above relationship. 
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D. Cardiac Output. There are a variety of uncertainties regarding 
the accuracy of measurement of cardiac output in humans. Ex- 
tensive data on the changes in cardiac output in a wide variety of 
physiological and pathological conditions was accumulated by A. 
Grollman (1932) using the acetylene method. More recent work by 
A. Cournand e¢ al. (1945) and L. Werko e¢ al. (1949a) using the di- 
rect Fick method with venous blood obtained by catheterization of 
the right heart gives values averaging 25 per cent higher than 
those obtained by Grollman. These authors believe their method 
to be free from some of the errors inherent in Grollman’s tech- 
nique. The dye dilution technique (review by Dow, 1956) gives 
cardiac output values close to those of the direct Fick as shown 
by Werko e¢ al. (1949b). It is the opinion of the present author that 
a precise evaluation of the accuracy of any of these methods in 
the hands of any particular group of investigators should be based 
on a comparison with direct measurement of cardiac output in 
heart-lung preparations. 

For this analysis, only cardiac outputs measured either by the 
direct Fick or dye dilution methods will be used. 


E. Arterial Blood Pressure. The conventional clinical method of 
determining blood pressure consists of inflating an air-filled cuff 
about the upper arm, then slowly deflating it and listening with a 
stethoscope distal to the cuff, noting the appearance of the sound 
which is taken to indicate systolic pressure and noting the muffling 
or disappearance of the sound which is taken to be diastolic pres- 
sure. In a study using 47 human subjects, brachial blood pres- 
sures were determined almost simultaneously by the clinical 
method and by arterial puncture using strain gauge and capacitance 
pressure transducers (Roberts e¢ al., 1953). For systolic pres- 
sures, values obtained by the clinical method averaged 12 mm Hg 
lower than those obtained by using the capacitance manometer and 
about 4 mm Hg lower than values obtained by using the strain- 
gauge manometer. For diastolic pressures, the clinical method 
using the muffle criterion gave values about 3-4 mm Hg higher 
than those obtained by direct arterial puncture using either type of 
pressure transducer. Diastolic pressure as determined by the dis- 
appearance of sound, as is sometimes done, gives values about 7 
mm Hg lower than those obtained by direct arterial puncture. 
Hence, it seems that the common clinical method of blood pressure 
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determination agrees reasonably well with direct measurement by 
arterial puncture. Unfortunately, in many published blood pressure 
studies, it is not stated whether muffle or disappearance of sound 
was used as the criterion of diastolic pressure. Other studies 
similar to the above suggest that muffle may be inferior to dis- 
appearance of sound as a criterion of diastolic pressure (Bordley 
et_al., 1951). 

Aortic pressures are nearly equal to those recorded in the radial 
artery (Wiggers, 1952). 


F. Venous Pressure. Venous pressure varies considerably in 
different parts of the body. The average venous pressure in the 
median basilic vein of the arm of the normal adult is given as 97 
mm of water (Burch, 1950), which equals 7 mm Hg, the value that 
will be used in this analysis. 


Mathematical Analysis of the Arterial System as a Whole 


Because the parabolic approximation of the aortic pressure- 
volume curve does not hold at pressures above 170 mm Hg, it is 
necessary to analyze separately Case I where P_ <170 mm Hg and 
Case II where P. > 170 mm Hg (P, = systolic arterial pressure), 


Case I: 1. Assumptions and Approximations. 
a. The mathematical model of the arterial system is as shown in 
Figure 3. 
b. Pressure in the elastic reservoir is everywhere equal at all 
times (i.e., pulsewaves down the aorta are ignored). 
c. The volume of the elastic reservoir at all times is given by 


V =aP? + bP +c (if P< 170 mm Hg) (3) 
or 


V =a (170)? +b (170) +c =V,, (if P>170 mm Hg), — (4) 


170 
where: V = total volume of the elastic reservoir, P = pressure in 
the reservoir. 

d. Flow (F) from the arterial reservoir to the venous reservoir 

per minute is given by the following equation: 
(P -P,) =RF (5) 

where: 

P, = venous pressure 


FR = constant, which is usually called total peripheral resistance 
p =the instantaneous arterial pressure 
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FIGURE 3. Model of arterial system used in this analysis. 


e. With each heart beat, when the aortic valve opens, the pres- 
sure and volume in the elastic reservoir each rise immediately 
to a certain maximum value, which is maintained constant un- 
til the valve closes, after which the blood stored in the 
reservoir is gradually discharged through the peripheral re- 
sistance allowing a gradual fall of aortic pressure and volume. 
When the aortic valve opens again for the next heartbeat, the 
same sequence of events is repeated. The maximum aortic 
pressure is called systolic pressure (P.), and the minimum 
pressure is called diastolic pressure (P,). The time when the 
aortic valve is open is called systole (duration = T.), and the 
time when it is closed is called diastole (duration = T,). The 
time course of aortic pressure as described above is shown in 
Figure 4. 

Actually the exact form of aortic pressure as a function of 
time during systole varies in different parts of the aorta. Pre- 
cise information on how this function varies in different situ- 
ations is not known. The straight-line approximation through- 
out systole was used to simplify the analysis. 


2. The Analysis. 
During diastole, from equation (5), 
aye* 4 


ete hae (6) 
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FIGURE 4. Approximation of the curve of aortic pressure vs. time 
which is used in this analysis for P < 170 mm Hg (Case I). 


By differentiation of equation (3) with respect to time, we also 
have 


ae = (2aP + 6) a (7) 
dt dt 
Combining equations (6) and (7), we have 
Ep Wye ey (8) 
dt R ; 


Separating variables and integrating over diastolic time (Fig. 


4), we find 
d (2aP +b) +b) T ota 
Jac carpe ates lag (9) 


Performing the integration, we have after simplification: 


P= he 

2a(P, — P,) + 2.308 (2aP, + 6) log,,| S—* | = 4. (10) 

: Py ay Ps 

The usefulness of equation (10) will be discussed later. We 

shall now derive another important equation by starting with the 

fact that the total flow per heartbeat is the sum of the flow through 

total peripheral resistance during systole plus the flow during 
diastole. We have 

s 


S = 


aay ck (11) 
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where 


S_, = flow per beat 
V_, = flow per minute (cardiac output) 
N = heartbeats per minute 


P_-P 
From equation (5), flow during systole = — R Sealy Jo (ko) 
Flow per beat in diastole = change in volume of aortic reservoir 


from beginning to end of diastole = ceo! ee 
where 


V, = volume in systole and at beginning of diastole 
V, = volume at end of diastole 


Since total flow per beat equals flow in systole plus flow in di- 
astole, 
2 {ft PY 


Se SL SG (18) 


Eliminating V, and V, by using equation (3), we have 


¥. Af gPi= P) pace 
ee ee) (14) 


Equations (10) and (14) give us two equations relating systolic 
(P,) and diastolic (P,) arterial pressures, cardiac output (V_), 
total peripheral resistance (R), and the constants of aortic elas- 
ticity (@ and 6). Duration of systole (T,) and of diastole (T,) as a 
function of heart rate can be obtained from Figure 2. Equations 
(10) and (14) can be solved simultaneously for the values of any 
two of these variables in terms of the other variables. 


Case II: 
If P, > 170 mm Hg, there must be an instantaneous drop in pres- 


dP 
sure to 170 mm Hg at the end of systole because pe (see Fig. 


1), and dV is finite because of continuing outflow. Hence, aortic 
pressure over time is as shown in Figure 5. 

The analysis of pressure and flow in diastole alone, as was 
done in the first part of the analysis of Case I, proceeds in exactly 
the same fashion, except that in equation (9) the lower limit of 
integration is 170 instead of P.. Proceeding the same as in Case 
I, we can then derive an equation corresponding to equation (10), 
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i 


FIGURE 5. Approximation of the curve of aortic pressure vs. time 
which is used in this analysis for P > 170 mm Hg (Case II). 


namely, 


. 51 Wt! ences oF 

2a (170 = P,) + 2.303 (2aP, + ) O8 10 i, Ps = ps: (15) 

The analysis of total flow throughout both parts of the heartbeat 

follows the same pattern as for Case I, except that the total flow 

in diastole equals change in volume from beginning to end of di- 

astole =V,,,-V,. Systolic flow is given by the same expression 

as in Case J. Continuing the analysis as in Case I, one derives 
the equation corresponding to equation (14), namely, 


| alta ioe i 

STERN Pai Aiiiak din! (AS es! oad BS 6 (18 a9 & (16) 
Equations (15) and (16) can be solved simultaneously for any pair 
of variables the same as can equations (10) and (14). 


Results 


A. Comparison of Theory with Experiment. Errors. M. Brandfon- 
brener e¢ al, (1955) and Landowne e¢ al. (1955) determined systolic 
and diastolic blood pressures and pulse rates with a strain-gauge 
manometer connected to the brachial artery and cardiac outputs 
(by dye dilution) almost simultaneously on 68 normal men of a 
variety of ages. Mean values were computed for different age 
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groups. Using this data and the mean aortic elasticity constants 
for excised human aortas calculated from the data of Remington et 
al. (1948), mean values for different age groups were substituted 
in equations (10) and (14). It should be noted that Landowne e¢ al. 
(1955) used values of FR calculated to be mean arterial pressure 
(P) divided by cardiac output, which assumes venous pressure to 
be equal to 0. Since our best estimate of mean venous pressure is 
7 mm Hg and since this value has been used in all computations 
in this paper, the measured values of R shown in Tables 2 to 5 are 
P=-7% 


, Where P and V_, are from the data in Brandfonbrener e¢ al. 
(1955) and Landowne e¢ al. (1955). Values of T, and 7, as a func- 
tion of N were obtained from Figure 2. The results of the compu- 
tations were compared with experiment. 

The first test made was to substitute V_, VN, &, a and 6 in equa- 
tions (12) and (14), solve for P, and P,, and compare with meas- 
ured values of P. and P, as shown in Table 1. As can be seen, 
the errors in P, are all less than 11 percent, and if the computed 
values of P. are compared with mean values of P, from arterial 
puncture, the errors are less than 6 per cent. It is ‘evident that if 
the computed values of P. are compared with the maximum meas- 
ured values of P., the errors are larger. Since the derivation of 
equations (10) and (14) involved an assumption that pressure was 
constant throughout systole, it is entirely reasonable that P, 


TABLE 1 


COMPUTATION OF P, AND P; For AGE GROUP MEANS WITH V_, V, AND 
R FROM BRANDFONBRENER ET AL. (1955) AND LANDOWNE £T AL. 
(1955) AND ELASTICITY DATA COMPUTED FROM REMINGTON 
ET AL. (1948) 


Errors as 
PALO 
Measured 


Age Group Computed 
(years) (mm Hg) 


Pi=7 


34 FREEMAN W. COPE 


(computed) should more closely approximate mean P, over systole 


(P.) than maximum systolic pressure (P, max.). The author has 
found by plotting P, vs. P, max. for the data of Landowne e¢ al. 
(1955) that P, =0. 85 le max. Hence, B; can be computed from 


P, max. with reasonably. small error if ok the latter is available. 

“The accuracy of equations (10) and (14) can be studied in an- 
other way by substituting Landowne’s values for Ps P,, N, and 
age group means for @ and 6 from the data of Remington et al. 
(1948), and computing R and V_, which are compared with meas- 
ured values of R and V_ in Table 2. As can be seen, four out of 
eight errors are 15 per cent or greater, which cannot be considered 
very satisfactory. 


TABLE 2 


COMPUTATION OF # AND a For DIFFERENT AGE GROUPS 


Age Group 
(years) 


Py: Pas and WN from Brandfonbrener et al. (1955) and Landowne ef al. 
cfosss and a and d as in Table I. Pots Fe is in liters per minute and 
mm Hg X min. 


Fis in 
liters 


When the present author plotted a vs. 6 for all of the cases of 
Remington et al. (1948), it was found that the points came very 
close to all being on a straight line, so that one may write a = 
(-3.06 + 0.9) x 107~* with only small error. This equation was sub- 
stituted in equations (10) and (14) so as to eliminate a. The re- 
sulting equations were solved simultaneously for R and 6, and 
then for V_ and 6. Values from the data of Brandfonbrener e¢ al. 
(1955) were substituted and the computed values of R, V_, and 6 
compared with the measured values from the same group of pa- 
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tients. Comparisons of computed with measured values of V_ and 
s m 
F are shown in Table 3. All errors are less than 7 per cent. 


TABLE 3 


COMPARISON OF COMPUTED AND MEASURED VALUES OF FP AND V. 
m 


Errors 


Age Group : ‘ as % of 
(years) Resistance (f) Cardiac Output (V_) een 
ured 


Computed Measured | Computed Measured 


Comparison of computed and measured values of resistance for different 
age groups from Brandfonbrener et al. (1955) and Landowne e¢ al. (1955). 
The following relationships were also used in the above computation: 
P,=P_, P=, 4 = (3.08 + 0.9) x 10%. 


The values of 6, when computed in these two different ways, are 
compared in Table 4. It is evident that in five out of seven cases, 
the agreement is 3 per cent or better. However, if one set of the 
computed values of 6 is compared with values of 6 determined from 
isolated aortas (Remington et al., 1948) as in Table 5, it is seen 
that in three out of four cases the errors lie between 20-25 per 
cent. Since the values of 5 computed in the two different ways 
show such good agreement (Table 4) but do not agree well with 
measurements on isolated aortas (Table 5), and since much more 
accurate results are obtained in computing V,, and R when 6 is an 
unknown than when values of 6 are used, it seems likely that the 
value of 6 in the intact human is significantly different from its 
value in the excised aorta. This error in 5 apparently accounts for 
a large part of the errors in Table 2 and some of them in Table 1. 
It would seem that if errors in the measurements of 5b are elimi- 
nated, then calculations using the equations developed in this 
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TABLE 4 


CALCULATION OF b IN TWO DIFFERENT WAYS AND COMPARISON 
WITH EACH OTHER 


Calculated Values of } 


Age Group Sa pee ee eae % Difference between 
(years) Ff and b Vi, and 6 Calculated Values 
Unknowns Unknowns 


Data from Brandfonbrener et al. (1955) and Landowne e¢ al. (1955). 


paper may be expected usually to give accuracy within 5 per cent, 
since this is what was found (Table 3) when errors due to 4 were 
prevented by making 6 an unknown. 


B. Applications. An exact solution of equations (10) and (14) 
for P, and P, cannot be accomplished, but a solution by iterative 
methods can be done to any desired degree of accuracy, although 
the computation involved is tedious. The author turned this work 
over to the Aeronautical Computer Laboratory at the Naval Air De- 
velopment Center. P. and P, were obtained by the Newton-Raphson 


TABLE 5 


COMPARISON OF CALCULATED VALUES OF } (6 AND 2 UNKNOWNS) IN 
LIVING HUMAN WITH VALUES OF } MEASURED IN ISOLATED AORTA 


Age Group 6 Calculated in 6 Measured in Difference as % of 
(years) Living Human Isolated Aorta Measured Value 


Data from Remington et al. (1948). 
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method of approximation (Scarborough, 1930), with the computation 
being done on an IBM 650 digital computer. An additional portion 
of the computer’s program stated that if P, >170 mm Hg, then P. 
and P, were to be recomputed using equations (15) and (16) of 
Case il. 

Since it now became possible to compute P_ and P, easily and 
quickly, the way became open for a quantitative study of the ef- 
fects on P, or P, of changing any of the other variables singly, or in 
combination. The results of such computations will be presented 
in a future communication from this laboratory. 

The equations developed in this paper allow the calculation of 
the elastic characteristics of the aorta in the intact human, An- 
other publication from this laboratory will describe these applica- 
tions in detail and will demonstrate that exercise and certain 
drugs can cause significant changes in the aortic elastic charac- 
teristics. 


Summary of Some Important Restrictions in the Application of 
this Theory 


The relation between the physical variables of the human sys- 
temic arterial system is described by equations (10) and (14) if 
P, £170 mm Hg and by equations (15) and (16) if P, > 170 mm Hg. 
Either pair of equations may be solved for any two variables in 
terms of the others. Solutions for P, or P, can only be obtained 
by approximation methods, although one can solve for the other 
variables by standard algebraic manipulations. 

The value of P, in these equations must be taken as mean pres- 
sure over time during systole as measured by strain-gauge ma- 
nometer in the brachial artery. By plotting a scattergram, it was 
found that mean systolic pressure = 0.85 maximum systolic pres- 
sure (by strain-gauge manometer), so the former can be calculated 
from the latter. Systolic pressure measured by cuff is slightly 
lower than the above as discussed in Section E. Cuff measure- 
ments of P, are also somewhat inaccurate (Sec. E). 


’ (P = mean arterial 


The correct value of R is by definition 

m 
pressure). If R is obtained from direct measurement of P and V_, 
the same estimate of P, must be used as was used in the other 


calculations. 
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The relationship a = (—3.05 + 0.9) x 10 * was found to hold be- 
tween the constants of aortic elasticity, so that the constant a 
may be eliminated from the pair of simultaneous equations. Hence, 
they may be solved for 6 and any other variable. 

Values of } in the living and dead aorta are significantly differ- 
ent. It will be shown in another publication that in the living 
person, 6 may vary significantly from time to time. 
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ANALYSIS OF TRACER EXPERIMENTS: IV. THE KINETICS OF 
GENERAL N COMPARTMENT SYSTEMS* 


H. E. HART** 
DEPARTMENT OF PHYSICS, COLLEGE OF THE CITY OF NEW YORK AND 
DIVISION OF NEOPLASTIC DISEASES, MONTEFIORE HOSPITAL, 
NEW YorK, N.Y. 


The methods of C. W. Sheppard and A. S. Householder (Jour App. 
Physcis, 22, 510-20, 1951), H. D. Landahl (Bull. Math, Biophysics, 16, 
151-54, 1954) and H. E. Hart (Bull. Math. Biophysics, 17, 87-94, 1955; 
ibid., 19, 61-72, 1957; 2b¢d,., 20, 281-87, 1958) are employed in studying 
the kinetics of general NV compartment systems. It is shown that the 
nature of the transfer processes occurring in fluid flow systems and the 
chemical processes occurring in quadratic systems and in catalyzed 
quadratic systems can in principle be completely determined for all 
polynomial dependencies. Systems involving three-body and higher 
order interactions can be completely solved, however, only if supplemen- 
tary information is available. 


Introduction. The development of general tracer and combined 
tracer perturbation methods has made possible analysis of the reg- 
ulatory mechanisms of fluid flow systems (Sheppard and House- 
holder, 1951; Landahl, 1954; Hart, 1955, 1957, 1958t). It is char- 
acteristic of fluid flow systems that direct interaction between 
constituents of different compartments cannot occur. Regulation 
in such a system must be essentially valvelike. The passage of a 
constituent from a physical compartment 1 into another physical 
compartment 3 (i.e., X, —+ X,) may be regulated by the same or a 
different compound in compartment 2 (i.e., X,), but X, as such 
does not combine with X, to form X,. This lack of interaction 
between constituents is necessary in order for the Qiig of equa- 
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(30-1)-1551. 
**Presently National Science Foundation Science Faculty Fellow, 
Department of Biophysics, Yale University. 
tHereinafter called references I, Il, and II, respectively. 
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tions (1) in reference II to be equal to the transfer functions a; 
of equations (5), reference Il. 

If the components of the system do interact, then the number of 
compartments may be thought of as given by N = PC, where P is 
the number of physical compartments and C is the number of chemi- 
cal substances involved in the system. V, =V,,V,,V,...Vp are 
the compartment volumes (assumed constant in this treatment), and 
X,=X,,X%,...Xy are the concentrations expressed in grams, 
moles or molecules per unit volume. For incompletely character- 
ized materials, such as certain enzymes, X, may be expressed as 
activity per unit volume. If not all substances occur in each physi- 
cal compartment, then some of the concentrations X, will always 
vanish. If a substance is localized within a compartment or ex- 
cretion occurs, no new compartments need be invoked, the processes 
being treated as in reference I. 

In terms of the above notation, it can be shown that tracer ex- 
periments carried out on systems with interacting components can 
give rise to ambiguities. For example, the formation of X, as il- 
lustrated in Figures la and 1b below can satisfy the same equation: 
X, = aX?, (1) 

In Figure 1a two molecules of the material X , combine to form 
X, and a is proportional to the probability of formation. In Figure 
1b, however, the term aX, corresponds to the type of processes 


X X\ X) 


aX, 


FIGURE la FIGURE 1b 
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treated in reference II]. Thus, the mechanism regulating the ‘‘flow” 
of X, into X, simply has the functional dependence aX ,. 

If tracer experiments with labeled X , are carried out on the sys- 
tem of Figure 1b, a resulting equation, as expected, will be: 


Ry, = (aX ,)X 4) (2) 


For the system of Figure la, however, there is twice as much 
X, in each molecule of X,, and this gives rise to the equation: 


OME CIES ae Oe (3) 


From the point of view of reference III, equation (8) would there- 
fore suggest a valvelike regulatory mechanism = 2 aX,. 

A similar inconsistency occurs in the systems illustrated in 
Figures 2a and 2b, where in each case the formation of X, is 
given by: 


X= aGX,4,. (4) 
X Xo Xj 
[nteraction Valve 
Fegulation regulation 
OXo 
Xz Xz 
FIGURE 2a FIGURE 2b 


Experiments in which X, and X, are labeled with tracers give 
rise to the result for Figure 2a:* 


X,, = (@X,)X1, + (aX ,)Xq,-** (5) 
For the system of Figure 2b: 
Kp (akg Xyy (6) 


*In some cases only a portion of a given molecule X, will be employed 
in actually forming X3. Therefore, depending upon the position of the 
radioisotope label, X; will exhibit valve or interaction regulation. 

**The form of equations (4) and (5) is preserved if X3 is a polymer of 
X1,Xq. Note that since tracer concentrations are by definition small, 
multiply-labeled individual molecules can be neglected. 
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In the present work, the use of tracer and perturbation tech- 
niques in the analysis of general systems involving valve-like 
regulation (i.e., enzymes, catalysts, and fluid flow processes) 
and interaction regulation will be discussed. 


Table of Notation 


AG (¢) = concentration of material in the jth compartment which is 
regulating in a valve-like manner the passage of material 
in the ith compartment; associated with this concentration 
will be an exponent «’. 

‘X, (¢) = concentration of material in compartment j which is flow- 
ing into another compartment 7; associated with this con- 
centration will be an exponent 4%. If «’*> 1 or if more 
than one ‘‘vector’’ concentration appears in a single term, 
an interaction process is involved.* 

X, (¢) = concentration of material in compartment ¢. 

X,, (¢) = concentration of freely diffusible material in the zth com- 
partment labeled with the /th species of isotope. 


I. General Theory 


It follows from the Table of Notation above and equations (5) of 
reference III that using perturbation methods alone one obtains:** 


tot De Oa aie ee 


z J 
oa > Aree SSE 
Note that in equations (7) all X’s appear symmetrically without 
preferentially specifying any particular X,.. This is appropriate 
when treating interactions between different chemical constituents. 
Otherwise, as an example, the simple process X, + X, — X, in- 
volving transfers from both compartments 1 and 2 into compartment 


*It should be noted that in the definitions of iX; and 'X; the pre-super- 
script t and the boldface are used to indicate that the concentration, Xi 
arises in a certain way in equations involving X;. The additional sym- 
bols do not change the value or meaning of X;, however. 

**For fluid flow systems: CiaB...w= Bita- 1,B...0+ Bi20, B-1...0 + 


BinaB...0-1 =e: » Bijak..e- 1...@ (Compare eq. 5, ref. I). 
i 
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3 might be counted ¢wice. Equations (7) of course do not con- 
pletely reflect the nature of the processes (i.e., whether valve- 
like or interaction) which are regulating the concentrations Ge 

In order to examine this latter problem, it is perhaps’ convenient 
to refer to the hypothetical reaction: 


Dae ¢ 
eee ee A) Fe Xi (8) 


4 

Equation (8) represents 2 ‘‘units’’ of X, and one unit of X, com- 
bining to form one “‘unit’’ of X,. The material X, in this example 
Serves not only as a constituent of X, but for generality is as- 
sumed in association with X, to catalyze the reaction as well. It 
follows then that the rate of formation of X, is determined by the 
product of three factors: the square of the concentration of X,, the 
concentration of X,, and the function F (X,,X,) which is assumed 
to characterize the catalytic effectiveness: 


aX , 
ms 


(9) 


If X, and X, are now labeled with an isotropic tracer, the con- 
centration of labeled X, arising from this particular reaction will 
be given by: 

aX 4 
dt 


EE PONG Xf: (10) 


where the three terms within the braces represent the three pos- 
sible means by which X, can be labeled. Since in general labeled 
elements are relatively rare, terms such as X,,X,,X, or X,X,,X 9, 
are omitted as in equation (5). 

Rewriting equation (10): 


aX 41 


{2X XX, + (4)? X13 (11) 


1? 
it is seen that the factor ‘‘2’’ arises from the fact that 2 units of 
X, are present in each molecule of X,, so that for the same spe- 

mite activity of X,, and X,, there is twice as much likelihood of 
X, being labeled in an x, eaten as by an X, element. Note, 
peravse that the dependence of the catalyzing Ainetion F upon X, 
X, does not affect the statistics of relative labeling Soa 
abilities. 


46 H. E. HART 


Equations (10) and (11) can be rewritten in a form somewhat 
more suitable for generalization: 


dX ,, F) J | 

= ea ee el Xs —— [(X,)°X,] 12 

dt Fu Xe) bear ax [( Ate ate ee ( p 2 ( ) 

a a x, ? F(X, X,) X? X,I, 
j=1,2 

where 
JX OX, 
te Oar = Bae (13) 
OX, OX, 


and the boldface serves to emphasize that an actual “‘flow”’ of X, 
or X, into X, rather than catalytic action is involved. Thus equa- 
tion (12) requires that only the interaction regulation factors be 
differentiated. Since the constituent X, may be formed in several 
ways with different initial components and different catalyzing 
functions it is perhaps appropriate to note explicitly that Xe is 
defined to be the concentration of a substance X. which in a given 
process is catalyzing the formation (or degradation) of X,; also 
xX, is defined to be the concentration of a substance which in a 
given process is being incorporated as a constituent of X,. 

If the further assumption is made that the catalyzing functions 
can be represented as simple polynomials in the concentrations 
X,, then it follews that for a general system: 


‘ 0 
Xin 7 »: Xi C 0? B?...00°; Ott Bi? eGlae d(?X.) 


i, BY...05 OO BP co! i 
ee re ee ee ee (14) 


where &’, B”...@° are the exponents arising from the polynomial 
expansions of F; ’*.8’’...@** are exponents reflecting the num- 
ber of units of each substance which are employed as constituents 
in forming X, as indicated in the equation below: 


Ot AE BR ORT A Mate a (15) 


and the constants C,,,, Bo veaco”s 07 B77... 92%, ate determined by the 
polynomial expansions of the catalyzing functions for the reaction 
of equation (15). 
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Carrying out the differentiations in equations (8): 


X + ) x 
HC eT sir, Sen CdS Hee ag 


L 7 
Xx) fa G 


Bete a Ose Os ioe ta ennai Ole cee 
Qa , _ 
Me eet el ee X=.0,. S16) 


where & = %” + &’” etc. Rewriting equations (16): 


a7 Br... £77...0/* = 0 


Coro (18) 


£5(& ~ eS lo hkk fies S— C72, erento > w’?; oR eeey SIA ee, (Nk prt p(B 


It follows from equations (16) or (17) that the functional de- 
pendence of the individual transfer parameters between compart- 
ments can therefore be determined by using the perturbation-tracer 
methods of reference II (i.e., CiiaB... ©... » can be calculated).* 
From equation (18) it can be seen, however, that the nature of the 
regulation giving rise to the observed functional dependence (i.e., 
whether valve or interaction regulation) will often remain unde- 
termined since the individual C; y= yor w—wet; 02%, BP"... wee CAN 
not always be evaluated. 


II. Special Systems 


1. Fluid Flow Systems. For fluid flow systems e€’’ = 1, all other 
doubly primed exponents vanish and equations (11) above reduce to 


Pe aN aX so XE Xyim 0, (19) 
Fai sie 


*A simplified method of evaluating coefficients, owing to Prof. H. D. 
Landahl, appears in the Appendix. 
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corresponding to the treatment in reference II, and the B., Py es 
can therefore be found. 

2. Quadratic Systems. For systems involving only fluid flow and 
two-body interactions: 


N N 
ket YS Ay Ak, > Yo Bi AAG =” (20) 
j=l i,k = 1 
where A,, and B, ., are coefficients which can in principle be de- 


termined by Soteeon methods alone. By invoking the mathe- 
matical formality of an additional stable concentration X, = 1, the 
linear terms need not appear explicitly and the appearance of equa- 
tion (20) is simplified: 


N 
X,+ Da Cape dp hgimicak ods (21) 
j,k = 


When X, is formed it arises presumably from processes involving 
either the incorporation of X. or X, or both. Indicating these three 
possibilities explicitly equation (21) becomes: 


N 


Ge: Sind Naa Levies), a? ak ot ree), ee 
i,k = 
where the subscripts of the coefficients C;, S[k? C,, + pk? C; -jk Serve to 


identify the processes with which they are associated. ‘Comparing 
the notations of equations (22) and (14): 


C. 1 See =C,, :000...0; 110...0 C i 12 o Ce 100...0;010...0 


(23) 


C i312 3 Cote 100...0 Cin = Ci ;000...0; 200...0° 


It follows from equation (22) that the equations for tracer studies 
become: 


N 
Rigs): El Are Gent Ry dae a 
j,k =0 


(Coie + Cazj,) X, Xj} = 0- (24) 
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The coefficients 
ik? - ef 
t3j i;/k i;jk (25) 
eet 
poke = Cistk + Ci spk 
can be determined numerically as indicated in the Appendix. Since 


Citi = C sik Bs Cistk 43 C sik (26) 


is known from equations (7), it follows that the coefficients Cire 
Ci ik and Ci ig which completely determine the nature of the regu- 
lation can be evaluated: 


i;[k i;jk bsgeul a 
(27) 
C iste a C sik 6 Ci ejK? 
C i sik a Ci ik fee Ci irk : 


3. Quadratic Enzyme Systems. Consider the system of equations: 


; iy okys i 
te le A AS (28) 
jkz 

where “Sa represents any particular product of a number of enzyme 
concentrations known to be catalyzing the conversion of X, or X, 
into X,. If regulatory terms are restricted to *X, Xecans x, ‘Ak, 
and *X °K, ‘X,, then C istkz? Cpe Cisikz can be evaluated using 
the methods outlined for simple quadratic systems: 


C istke 7 Be pcalen + Ciieke = Ciik ee 
Ci ike ~ Ci sjkz = C isin? sz (29) 
Ores 7S ssyee ~ Cash e” 


4. Cubic Systems. Third-order systems of the form: 


Xt Dy Cram X% Xf IXY = 0; la] + 181 + ly! =3, 80) 
a By jkm 
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where all constituents can participate directly in flow processes 
(i.e., ‘Xi, xe UXea # 0) cannot in general be solved unless ad- 
ditional information is available. Expanding a typical term (i.e., 


= Bi=y =1): 


Ky Oc XOX eT Oe ee ae 

#0 pack) Xp Xe Ue ee eee (31) 
sire ee 7, ee Ges Ep GRD 

POR ae hy A ae 


Since only four parameters C, Sih C.. jkm? Ceikem and C, é 
are determined in the method outlined,’ oe follows that comiiete 
characterization of third- and higher-order systems is not possible 
at the present time without making use of supplementary informa- 
tion. The incompleteness of the analysis of third- and higher-order 
systems would seem to reflect the fact that separate measurements 
of the possibly intermediate two-body interactions are not being 
employed. It is as if these intermediate compounds constituted 
inaccessible compartments. Higher-order systems are therefore 
often incompletely solvable not because they are inherently dif- 
ferent from quadratic systems but simply because measurements 
adequate to reduce them to the form of two-body or two-body-plus- 
enzyme systems are not available. 


III, Discussion 


From the foregoing it is apparent that in principle and under 
ideal circumstances the kinetics of most biochemical systems can 
be determined. In actual practice, however, the experimental ac- 
curacy possible is often far below that required for effective use 
of the formalism. Intelligent use of any of the methods outlined 
therefore requires that supplementary information be employed in 
such a way as to cut down on the variables and minimize the ef- 
fect of errors in measurement. 

It would seem that equations (7) bear some resemblance to the 
equations arising in the study of dynamical systems. Whether 
these analogies can be developed into a Lagrangian or Hamiltonian 
type treatment of tracer experiments in multi-compartment systems 
and, if so, whether any essentially new results can arise from such 
considerations remains to be seen. 
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their interest and support of this work. I should like to thank Dr. 
Robert Hatcher for his careful review of the theory and calcula- 
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APPENDIX 


Assume that it is known a priori that the summations in equa- 
tions (7) contain only a few terms and also that simultaneous meas- 
urements can be carried out in all compartments. Then equations 
(7) can be written explicitly at many different times 2: 


Xt) Cia B.ww Xs LENG = 0, (32) 
aB...@ 
and the coefficients C, __» evaluated. Similarly, rewriting equa- 


tions (24) explicitly for Mitoront times ¢: 


Xi, + a ree s Csi Area, + 
j,k 


[c 1X,,%j.,5=0. (88) 


i;jk C;, jk 
Carrying out measurements again at a sufficient number of times 
[C,. ie C sie and [Ci + C;,;,] can also be determined. 

Only a single isotope is required to carry out the analysis of 
the system represented by equations (21) because the functional 
dependence is actually known to within a finite number of coef- 
ficients. The ‘‘amount’’ of information being obtained is therefore 
quite similar to that sought for steady state systems and the methods 
analogous [see equations (4), (5), and (6) of reference I]. From a 
slightly different point of view, it is seen that only one isotope is 
required to analyze such an N compartment non-steady state sys- 
tem instead of the N-1 isotopes usually required because actually 
an infinite amount of information is already known—namely, that 
all but at most a given finite set of coefficients in Taylor series 


of equations (7) vanish. 
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Stress systems in a deformed medium are in general of two kinds: 
(i) dilation-producing and (ii) nondilation-producing. Ordinarily the dila- 
tion is computed from a knowledge of the stresses, which in turn are 
obtained by solving the general elastostatic equations of equilibrium, 
under specified boundary conditions. If, however, interest is centered 
only in dilation-producing stresses, or in the event that only such stresses 
exist in the deformed system, and the density changes in the medium may 
be experimentally determined or postulated a priori as functions of the 
coordinates, these stresses may be directly obtained by a simplified 
procedure of solving for a particular integral of the Poisson equation 
without appeal to boundary conditions.** In irrotational systems in which 
conservative body forces exist, there is a direct relationship between 
the dilation and the distributed body forces. Reference to both purely 
physical and biophysical systems is made. 


It is usual to compute the dilation, 6 = oe €;;, in a deformed 

i 
elastic medium from a knowledge of the stresses, which have been ob- 
tained by solution of the general elastostatic equation of equilibrium 


uV2u4+ (A+ p)V O(r) = - F(r). (1) 


However to solve these equations exactly, boundary conditions 
are required and their integration is often difficult. The displace- 
ment vector obtained, from which stresses and strains may be 
directly calculated through simple coordinate differentiation, is by 


*This investigation was supported by a PHS research grant, RG-5067, 
from the Public Health Service. 
**The solutions are unique to the extent that non-dilational stresses 


may be superimposed. 
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the Helmholtz theorem, the sum of a gradient and curl term, 
u=Vd+VxA. (2) 


The curl term is dilationless, since div curl = 0, and reflects only 
the displacements arising in the medium to satisfy the applied 
boundary conditions (Morse and Feshbach, 1953). The gradient 
term on the other hand, contains the internal displacements, 
stresses and strains, salely due to compression and expansion of 
volume elements within the medium. These may arise in part from 
applied surface tractions and in part from body forces such as 
gravity. 

There are problems in both physical and biophysical systems in 
which only dilation-producing stresses in the deformed medium are 
encountered. One can also visualize systems in which, although 
both dilation-producing and nondilation-producing stresses occur, 
interest is centered on the former. Thus an electric current gener- 
ating heat at a uniform rate in a straight solid conductor with no 
external restraints develops a thermal stress field of only the 
dilation-producing type (Wang, 1953). Similarly, a blood vessel in 
a state of contraction behaves like a rigid tube with respect to the 
uniform flow of fluid through it (Folkow and Birger, 1956). The 
pressure of the fluid, although a boundary condition itself, sets up 
a pure dilation-producing stress system. 

If now the dilation 6 in a deformed medium instead of being com- 
puted from the stresses is ascertained experimentally by measuring 
the density of the medium before and after deformation, or postu- 
lated a priori, one can carry out the reverse procedure, that of 
obtaining the internal stresses and strains which are dilation- 
producing. The recent development of rapid freezing techniques 
with helium, preserving molecular order, and of photometric density 
measurements, puts the reality of obtaining the density as a func- 
tion of the coordinates within tangible reach. One can then carry 
out the following analysis: 

Consider an element of mass, Am, in the medium. Because 
mass is conserved, we have 


Am = p(r)Av = p(r) Av’, (3) 


where p(r) and p’(r) are the densities and Av and Av’ the volume 
elements before and after deformation, respectively. Hence one 
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obtains for the dilation, (Sommerfeld, 1950): 


dye = (4) 
v p (r) 

Equation (4) expresses the dilation as a function of the coordi- 
nates of the medium in terms of the densities before and -after 
deformation, which in turn are also functions of the coordinates. 

Though the ascertainment of the dilation may be of interest in 
itself, the question arises whether any information concerning the 
body forces, stresses and displacements in an elastic medium can 
be derived from an experimentally secured (or a priori postulated) 
dilation. In other words: /f the dilation is given, what can one 
say about the state of an elastic medium? 

It will be shown first that, under certain circumstances, one can 
relate the body forces acting in an elastic medium directly to the 
dilation. 

We restrict ourselves to the consideration of homogeneous, 
isotropic media, for which the elastostatic equation of equilibrium 
is equation (1):* 


pV2u + (A+ yu) V A(t) = - F(2), 


where u is the displacement vector, 6(r) the dilation as a function 
of the coordinates, F(r) the body force per unit volume as a func- 
tion of the coordinates, and A and yp are the elastic moduli (Lamé 
Constants) of the medium. 

Substitution of equation (2) into equation (1) yields another form 
of the equilibrium condition, 


(A + 2n)V 0+ pV x (V2A) =—F(r). (5) 


From equation (5) it then appears that in principle the body force 
and the dilation may be specified independently of each other. 
There is, however, an important exception. If the body force is 
conservative, i.e., F =-— VV, where V is a scalar potential, equa- 
tion (5) can be rewritten as follows: 


V((A+ 2n)0-V]=8 (6) 
uV x (V72A)=- 8B, 


*Equation (1) is sometimes associated with the name of Navier Cf. 
I. S. Sokolnikoff, 1956, p» 73. 
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where the vector B, being at once the gradient of a scalar and the 
curl of a vector, is both solenoidal and lamellar. If it is known or 
assumed that the displacement u is derivable from a scalar poten- 
tial only, then from equation (2), Vx A=0. In this case, as well 
as in the case where the somewhat weaker assumption V x (V? A) = 
V2(V x A) = 0 applies, we have B=0. Hence, 


VU(A + 2n)9 - Vi = O, ie., (A + 2u) VO = VV =— F(r) or 
(A + 2n) O(r) = V(r) + constant. 


In other words, if V x (V2A)=0 and the body force is conserva- 
tive, equation (5) can be split into two parts: a ‘longitudinal’ part, 
i.e., the gradient part which involves the body forces and the 
dilation, and a ‘transverse’ part, i.e., the curl term. Both parts of 
equation (6) can then separately be equated to zero (cf. Morse and 
Feshbach, 1953, Vol. II, p. 1786). From equation (7) it can be 
seen that from a given dilation one can obtain the body force dis- 
tribution, and the converse: from a given body force potential one 
can obtain the dilation up to a constant. 

As an example one may consider the thermal stresses which are 
known to arise from a conservative body force in an elastic medium. 
Indeed, if there are no other body forces present, one can easily 
obtain the thermal stresses from the dilation, which in turn can be 
obtained by measuring the density of the medium as a function of 
the coordinates before and after heating, as was suggested above. 
The body force, from which the thermal stresses are derived, is 
given by (cf. Sokolnikoff, 1956, Sec. 99): 


F =—VV =-(8A + 2n) aVT(r), 


from which, the scalar potential is, by integration, 
V = (8A + 2u) OT + constant, (8) 


where & is the coefficient of linear expansion, and T7(r) is the 
temperature as a function of the coordinates. 
With the aid of equation (7) one then obtains the following rela- 


tion between the dilation and the temperature distribution in the 
medium: 


A + 2u 


rer 241) & 


- 0(r) + constant. (9) 
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It is thus seen that, when certain simplifying assumptions can 
be made, knowledge of the relationship between the distributed 
body forces and the dilation may be obtained. The relation be- 
tween the dilation and the elastic stresses and displacements is 
not so definite. This may already be inferred from the fact that 
stresses and displacements may exist in a deformed elastic medium 
without there being any dilation at all. The stress system usually 
designated as pure shear, for instance, is associated with a zero 
dilation. Thus, although one can compute the dilation directly 
from the stresses and/or displacements at the boundary, by Betti’s 
Method of integration (Sokolnikoff, 195€) the converse cannot be 
achieved; i.e., one cannot arrive at a unique stress system on the 
basis of just the dilation alone without further data such as 
boundary conditions. From the dilation only one can, however, 
arrive at a system of displacements and stresses, which is unique 
up to dilationless stresses and displacements. This procedure, 
besides being very simple, may be of value in situations where 
one is mainly interested in the dilation-producing part of the 
stresses and displacements, the nature of the dilationless stresses 
being of no importance.* 

It is convenient then to divide the elastic stresses into two 
types: dilation-producing and dilationless stresses. Similarly, 
deformations consist of two types: those associated with a non- 
zero dilation, and dilationless deformations. We propose to solve 
only for the dilation-producing stresses and deformations (dis- 
placements), which will be unique to the extent that dilationless 
components may coexist. The procedure outlined here amounts to 
finding a particular solution of a Poisson equation, or inquiring 
about the nature of a vector field, with only the sources, in this 
case the dilation, specified. The solution of such a problem is 
evidently non-unique, within the above indicated limits. 

Consider then the equation 


V.u= Ar), (10) 


*In this connection cf. Morse and Feshbach, loc. cit., p. 1786, where 
the following statement is made: ‘‘...the transverse parts of the solu- 
tion, M and N, are not particularly important in the static elastic case 
compared to the longitudinal part; they represent the unavoidable twist- 
ing of the medium in order to have the displacements required by the 
boundary conditions, but it is the longitudinal part which plays the major 
role with respect to internal strains and body forces, such as gravity.”’ 
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with the right-hand side specified. Using equation (2), equation 
(10) can be rewritten as 


V7 =r), (11) 
which is a Poisson equation, one particular solution of which is 
given by 

1 O(r’) 
$(r) = -— / See (12) 
4 volume |r oe r| 


In general one has to add a divergenceless solution 


U’= Vg'+VxA where V2g’=0 (13) 


to the displacements obtained from equation (12) by taking the 
eradient. This displacement is a solution of Laplace’s equation 
and is just the part necessary to fit the boundary conditions. 
Since, however, displacement (13) is dilationless, it is left un- 
determined in the procedure outlined here. It should also be 
noted that any two particular integrals of equation (11) can differ 
only by a harmonic potential, i.e., dilationless displacement. For, 
if #, 18 an integral of equation (11), and ¢, is another integral, 
V°5,-V°d,=V71d,-¢,]=V?S=0. From the displacements 
obtained above, the ‘significant’ stresses may be obtained by well- 
known procedures of linear elasticity. ‘Significant’ here means 
any system of stresses and displacements leading to the given 
non-zero dilation. 

The integral in equation (12) is to be taken over the volume of 
the medium in consideration, throughout which the dilation must 
be specified. It turns out, however, that once the integration has 
been carried out, the terms containing the surface parameters are 
harmonic and can thus for simplicity be rejected. This is indeed 
not surprising, since we do not expect the particular integral to 
depend on the boundary conditions. 

Consider, for example, a spherically symmetric dilation, varying 
exponentially in a spherical shell of inner and outer radii @ and 8, 
respectively. 


O(r) = O(r) = e”', (14) 


then, 
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The surface parameters a and 6b appear in the constants of the 
expression in the square bracket, but this expression is indeed 
harmonic (except at the origin), and thus gives rise to a divergence- 
less displacement. It may therefore be discarded in our con- 
sideration and the significant potential taken to be 


ey? 


Oe Pesce 2)- (16) 

This potential is entirely independent of the shape of the medium, 
and it will appear in any problem having the dilation specified in 
equation (14), irrespective of the shape of the deformed medium, 
whether spherical, ellipsoidal, etc. Any additional terms in the 
displacement, besides the one arising from equation (16), will be 


*Cf. Peirce, 1929, ##402, 403. 
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dilationless and necessary to meet any given boundary speci- 
fications. 

Since a solution of equation (10) is not necessarily a solution 
of equation (1) or (5), the question arises whether the displace- 
ment arising by taking the gradient of a particular solution of 
equation (11) can be made to satisfy the elastostatic equation of 
equilibrium. That this is so can be seen from the following con- 
sideration. For a specified 6(r), equation (11) yields a particular 
integral—which is a scalar potential ¢. By taking the gradient 
of this potential and comparing with equation (2), we see that any 
solution of equation (10) [or equation (11)] differs from a solution 
of Navier’s equation by at most a dilationless displacement. The 
curl term, being solenoidal, is left undetermined in this procedure 
as previously discussed, and therefore the statement that the 
equilibrium equation is, or rather can in principle be, satisfied is 
correct, although not very significant here. 


Summary 


Measuring or assuming @ prior’ the densities in an elastic 
medium before and after deformation, one can ascertain the dila- 
tion in the medium as a function of the coordinates. If one is 
justified in assuming that the curl part of the displacement (Helm- 
holtz equation) satisfies V?(V x A) = 0, then the dilation and body 
forces, which are assumed conservative, are directly related. If 
the primary interest is in the dilation-producing stresses and dis- 
placements in an elastic medium, these can be obtained directly 
by finding a particular integral of the Poisson equation (11), 
instead of solving the complicated elastostatic equations of 
equilibrium, subject to boundary conditions. Any such system of 
displacements and or stresses arrived at by the method outlined 
is, although obtained by a simplified procedure, not unique, that 
is, unique only up to dilationless components. Such cases may 
arise in either purely physical or biophysical systems. Recent 
technical developments may aid in obtaining density measurements 
as functions of the coordinates. 
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It is shown that the linear equations of the formal theory of mass be- 
havior proposed by the author in 1939 are first-order approximations of 
the biophysical theory proposed later. The second approximation con- 
tains cubic terms. It permits the expression of the time course of the 
change of behavior in terms of a quadrature. 


In previous publications (Rashevsky, 1939, 1947, 1959) we con- 
sidered a society composed of N individuals and divided into three 
classes, I, II, and III. The individuals of classes I and II, desig- 
nated as ‘‘actives,’’ determine their behavior because of inherent 
preferences. The individuals of class III, the ‘‘passives,’’ choose 
a given behavior ‘‘because most people do so.’’ Individuals of 
class I exhibit a behavior F,, those of class II a behavior Ry. 
The two behaviors are mutually exclusive. 

Let the number cf individuals in class I be X,, the number in 
class II, Y,. Then the number of passives is 


aes Ee Gey aa (1) 


Each class of actives influences the passives to accept the be- 
havior of that class. Denoting by X the number of passives which 
exhibit at any given time behavior F,; by Y the number of passives 
which exhibit behavior #,; and by a,, c,, and a, three constants, 
we assumed that the effect of the actives on the passives is de- 
scribed by the following equations: 


dX 

ie = 4X. + 4,X - Cg Y, - 4,Y; (2) 
dY 

ae = aX, - a,Y - CoX 6 -a,X. (3) 
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Because 
X+y=M, (4) 


we actually have only one independent equation. 
We have shown that when 


a,+a a, +4 ¥0 (6) 
) 1 0 1 

then X, regardless of initial conditions, becomes within a finite 

time equal to N’. In other words, when (5) holds, all passives ex- 

hibit behavior #,. Similarly, when 

Poss 


0 


N+ 


ao +a, ayo + a, 


Xo, (6) 


then all passives exhibit behavior &,. When neither (5) nor (6) is 
satisfied, then all the passives exhibit either behavior R, or Ro. 
If initially, (5) is satisfied and then at a given moment, due to 
changes in either X, or Y, or in a, and cy, Y, becomes satisfied, 
then a relatively sudden transition of the behavior of the passives 
from FP, to R, takes place. 

We also discussed the case where the constants a, and ¢), the 
so-called coefficients of influence, decrease with X, correspond- 


ingly with Y, so that, with a%, cx, €, and €’ as constants, 


ay = ak (1 - €X); Cy = cf(1—- & Y). (T) 
In that case equation (2) becomes 
dX 7 
se (2a, -agex, -c§fe&Y)X- 


[(a, — c* & Y)N’ —-akX,+c%¥]. (8) 


In this case the ratio X/Y varies continuously with X,/Y,, but 
only within a certain range. If X,/Y, exceeds a threshold value, 
X becomes equal to N’, X/Y =. 

Subsequently (Rashevsky, 1951, 1959), in attempting to give a 
biophysical interpretation to this type of behavior change, we con- 
sidered a population of N individuals, in which the tendency of 
each individual toward behavior R, or R, is distributed continu- 


ously. This tendency was considered as expressed by the differ- 
ence 


p=&,-& (9) 
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between certain central endogenic excitations €, and &,, where €, 
results in behavior F,, while €, results in behavior R,. If 6 = 0, 
the individual has an equal tendency toward Rk, and F,, and ex- 
hibits half the time behavior R,, half the time R,. If 6>0, the 
tendency toward FR, is greater, and the relative frequency of oc- 
currence of R, and R, is given by H. D. Landahl’s expression 
(1938). A corresponding situation holds for ¢ < 0. 
We assumed a continuous distribution N(d) of d, such that 


[ N@ae-n, (10) 


and considered that the effect on a given individual of seeing 
another individual exhibit behavior F, acts as a stimulus which in- 
creases €, according to the basic equations of the mathematical 
biology of the central nervous system (Rashevsky, 1960). We 
furthermore considered the simplest case, where this ‘‘imitative’’ 
effect is the same for all individuals. 

Denoting now by X and by Y the number of individuals which at 
any given time exhibit behavior 2, and R,, respectively, and by A 
and a two constants, we arrive at the following equation: 

d¢ 
ap 7 AA - N- 49. (11) 

If N(d) is symmetric with respect to ¢ = 0, having a maximum at 
that point, then (Rashevsky, 1959) X-Y, when plotted as a function 
of ¢, has the shape shown in Figure 1. Equilibria solutions of 
equation (11) are determined by the intersection of the X-Y curve 
with the straight line a¢. The slope of X — Y at 4 is propor- 
tional to N, and may be written BN. If a< ABN, then there are 


FIGURE 1 
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three equilibria, two stable and one unstable. For a > A3N, there 
is only one (stable) equilibrium at 6 = 0, X - Y =0. 

If N(#) preserves its shape but is displaced as a whole in the 
direction of negative 4, thus introducing a bias in favor of F,, 
then the X — Y curve is shifted to the right. If the society is at 
first [for symmetric N(¢)] in equilibrium 0,, then the majority ex- 
hibits behavior R,. As the bias toward 7, increases and this 
curve shifts to the right, a moment will come when the straight 
line will become tangent to the X — Y curve. At that point there 
will be only one stable equilibrium, 0,, and the majority will rela- 
tively suddenly change their behavior from f&, to R,. Thus we 
have here an analogy with the situation described by equation (8), 
namely within a certain range X — Y; therefore, X and Y vary con- 
tinuously, but at a certain point a sudden change occurs. 

In attempting to interpret equations (2), (3), and (8) in biophysi- 
cal terms, we suggested that they may represent approximations of 
the theory expressed by equation (11) (Rashevsky, 1959). The di- 
vision into three classes could be considered as a crude approxi- 
mation of a continuous distribution by a ‘‘rectangular’’ one. It 
was suggested that equation (8) may be a sort of second approxi- 
mation to the theory expressed by equation (11). The arguments, 
however, are not convincing. 

It is the purpose of the present paper to show that equation (2), 
for the particular case X, = Y, =0, is in an exact mathematical 
sense a first-order approximation of equation (11). In other words, 
the first-order approximation (11) leads to 


—=a,(X - Y), (12) 


which is just equation (2) for Xo = Y,=0. We shall also show 
that the next order approximation of the solution of equation (11) 
involves cubic terms and therefore that equation (8) can in no way 
be considered as a mathematical approximation to equation (11); 
although equations (8) and (11) describe two somewhat similar 
sociological phenomena. 

The exact shape of the curve in Figure 1 depends on N(¢). If 
the latter is a normal distribution, then the curve is expressed in 
terms of the integral of the normal distribution (Landau, 1950). 
Other assumptions about N(¢é) lead to different expressions 
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(Rashevsky, 1959). The general shape of the curve is, however, 
all that matters for our argument. 

If X - Y = F(d), where F has the shape shown in Figure 1, then 
the inverse function ¢ = F(X — Y) has the shape shown in Figure 
2. The slope of F at ¢ =0 is always finite. Hence the power 
expansion of F(X — Y) contains only odd powers. It also does not 


FIGURE 2 


contain any constant terms, if N(d) has a maximum at ¢ =0 (no 
bias). Thus the expansion has the form: 


pm E(k) aby (XIU) OX = Y) Ps es (18) 
But, since the slope of F at 6 = 0 is BN, 
1 
Seg nps 14) 


Thus, preserving only linear terms, we have 
X-—Y 
? = BN 


Por Key = Ne. (15) 


Introducing this into equation (11) we obtain 
d(X - Y) 


Fp (ABN - a)(X - ¥). (16) 
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Remembering that X + Y = N, and putting 


ABN-a 
a; (17) 
we obtain equation (12), Q.E.D. 
Equation (12) may be written 
a ee ee. (18) 


If a,>0, then for X <N/2 we have dX/dt<0, and the state 
X =0, Y,=N will be established. If X > N/2, then dX/dt>0, 
and the state Y =0, X = WN will be established. Thus for a, >0 
there are two possible configurations. If, however, @,< 0, then 
the only possible configuration is X - Y = 0. But the requirement 
a,>0 implies, from equation (17), @< ABN. As we have seen, 
this is the condition for the existence of two stable equilibria in 
the theory expressed by equation (11). The inequality a> ABN in 
that theory implies only one stable equilibrium for X -— ¥Y =0. 
Thus this feature is completely preserved in the first-order ap- 
proximation (12). 

The reason that as a first-order approximation to equation (11) 
we obtain equation (12) instead of the more general equation (2) 
lies in the circumstance that in the theory expressed by equation 
(11) there are no absolutely active individuals. An active indi- 
vidual in the sense of equation (2) is one with either d =+0 or 
@=-c. Their behavior never changes. Sut in the theory ex- 
pressed by equation (11), any individual with a finite ¢ is affected 
by imitation, 

In order to obtain equation (2) as a first approximation from the 
biophysical theory, we must introduce explicitly into the biophysi- 
cal situation completely active individuals. We may consider Xo 
individuals with ¢ = +0 and Y, individuals with 6 =-«. Thus 
we shall have two completely active classes and one class of in- 
dividuals in which ‘‘passivity’’ is continuously distributed. We 
Should now put, instead of equation (10): 


+00 
[ N(p)db =N’ =N-X,-Y). (19) 
Instead of equation (11) we now have 


dd 
Fp TAK - V) + aX q - eg Vy - a6. (ee 
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It is not surprising that the constant term 4X9 -— CoY, in equa- 
tion (20) plays a similar role as the variable term X - Y, in spite 
of the fact that the latter represents pure imitation, the former an 
active influence. Whether we do a certain thing because we pas- 
sively imitate someone else, or whether we do it because we are 
in some way actively influenced, in both cases the central excita- 
tion which causes the performance of the given act increases. 
Hence ¢ increases. 

Using again equations (15) and (17), we now obtain equation (2) 
as the linear approximation to equation (20). 

Now let us consider the second approximation to equation (11) 


and preserve in expansion (13) both the linear and the cubic terms. 
We now have: 


X-—Y 


to ey = 3 
Piste wa (Xi Y)” (21) 


Introducing this into equation (11) we find: 
d(X — Y) 
dt . 
(ABN - a)(X - Y) - ab, BNiA— Y)o (22) 


[1+ 35, BN(X - Y)?] 


The expression in brackets is always positive. Hence d(X - 
Y)/dt has the same sign as the right side of equation (22). 

If a< ABN, then the right side of equation (22) when plotted 
against X — Y is represented by the curve of Figure 3. It is seen 
at once that equation (22) admits three equilibria solutions, of 
which two (0, and 0,) are stable and one (0) unstable. We have 
the same situation as in the exact case represented in Figure 1. 
The approximate expression (22) has, however, the advantage that 
it can be directly integrated by separation of variables. 


KY 


FIGURE 3 
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If a> ABN, then the right side of equation (22) is always nega- 
tive for X — Y>0O and positive for X- Y<0. Hence again in 
agreement with the exact case, there is only one stable equilibrium 
at X= Y = 0. 

Equation (22) may be used to calculate approximately the neces- 
sary bias of ¢ in order to bring the society, say, from behavior Ff, 
to behavior R,. A shift of the curve of Figure 1 to the right corre- 
sponds to a shift of the curve of Figure 3 downward. When that 
shift is so large that the curve becomes tangent to the X - Y axis 
at the right maximum, then 0, will become the only stable equilib- 
rium. When the curve of Figure 3 is shifted downward, it becomes 
represented by the equation obtained by adding a constant 0, to 
the right side of equation (22). This gives, denoting X ~ Y by 2 
and the left side of (22) by u: 


u= 6, + (ABN - a)2 - ab,BN2?. (23) 


The abscissae of the maxima are given by du/dz = 0 or by 


eee Ae Nia ee (24) 
V 3ab,BN 


Since we are considering the right maxima, we use the positive 
sign. With this value of z, uw must be zero if the curve is to be 
brought to the z-axis. Hence, introducing equation (24) (with the 
positive sign) into equation (23) and equating the result to zero, 


we find: 
29 |(ABN -a)® 
b= = Seen 2 4 
Ota Ves ab,BN ao 


This is the value of the bias looked for. The time course of 
X — Y from 0, to 0, can now be obtained by integrating equation 
(22) with the initial condition that X - Y is given by equation (24). 

It may be remarked that in equation (20) the equilibria states are 
now determined by the intersection of the X — Y curve with the 
line 2=a¢ - (aX, - Cy)Y,). This line does not pass through the 
origin, but intersects the ¢-axis at the point (@)X 9 - Co Y,)/a. 
Hence the term a,X, - c,Y, has the effect of shifting the ¥ — Y 
curve along the ¢-axis; in other words, it has the same effect as a 
shift of N(d). 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A, Abbott Memorial Fund of the University of Chicago. 
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The principle of biotopological mapping (Rashevsky, 1954, Bull. Math, 
Biophysics, 16, 317-48) is given a generalized formulation, as the prin- 
ciple of relational epimorphism in biology. The connection between this 
principle and Robert Rosen’s representation of organisms by means of 
categories (1958, Bull. Math. Biophysics, 20, 317-41) is studied. Rosen’s 
theory of (9,8)-systems (1958, Bull. Math. Biophysics, 20, 245-60) is 
generalized by dropping the assumption that only terminal M; components 
are sending inputs into the #; components. It is shown that, if the pri- 
mordial organism is an (¥%,8)-system, then the higher organisms, ob- 
tained by a construction well discussed previously (1958, Bull. Math. 
Biophysics, 20, 71-93), are also (9,R)-systems. Several theorems about 
such derived (9,8 )-systems are demonstrated. 

It is shown that Rosen’s concept of an organism as a set of mappings 
throws light on phenomena of synesthesia and also leads to the conclu- 
sion that Gestalt phenomena must occur not only in the fields of visual 
and auditory perception but in perceptions of any modality. 


In a previous paper (Rashevsky, 1954) we proposed to encom- 
pass the mathematical description of important relational aspects 
of biology by representing an organism either by a topological 
complex or by a topological space. The basic feature of relational 
similarities between organisms was then formulated as the prin- 
ciple of biotopological mapping, Numerous applications of this 
principle were discussed in subsequent papers (Rashevsky, 1955; 
1956a,b). In the first paper (Rashevsky, 1954) we chose to repre- 
sent an organism by a directed graph which represents the physio- 
logical organization chart of the organism. Subsequently (Rashev- 
sky, 1958) we used discrete topological spaces. Still later (Ra- 
shevsky, 1959a) we found that we can deal directly with sets, 
which were not spaces, but which were characterized by certain 
relations between their elements. However, when both topological 
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spaces and sets are used, there always appears a directed graph 
associated with the corresponding space or set. Thus the theory 
of graphs seems to provide a useful tool for the description of the 
organism. 

In our original approach the vertices of the graph represented 
what we called ‘‘biological functions’ or ‘‘biological properties’’ 
of the organism. They are exemplified e.g. by locomotion, sensi- 
tivity, digestion, etc. The directed lines (arrows) of the graph 
stood merely for the binary relations: ‘‘The manifestation of a 
biological property A is directly followed by a manifestation of 
biological property B.”’ 

Instead of having the vertices of the graph represent the ‘‘bio- 
logical functions’’ or properties they may just as well represent 
the material, anatomical structures which are the vehicles of ‘‘bio- 
logical functions’’ or properties. The arrow in that case indicates 
that the functioning of a given structure (organ) results directly 
in a functioning of another structure. 

The inadequacy of the simple graph-theoretical representation, 
in spite of some of its successes, became apparent at once, Thus 
if a point A receives two arrows, one from a point 8, the other 
from a point C, it is not obvious whether A will function only if 
both B and C function, or whether either B or C are sufficient 
(Rashevsky, 1954, p. 338). This ambiguity was discussed in sub- 
sequent papers (Rashevsky, 1956a,b). Some artifices were intro- 
duced to solve the difficulty in certain cases (Rashevsky, 1956a), 
but no general way to avoid the ambiguity was found. 

Robert Rosen (1958a, hereafter referred to as RI) gave a some- 
what different interpretation to the graph-theoretical representation 
of an organism. Subsequently Rosen (1958b, 1959; hereafter re- 
ferred to as RII and RIII) succeeded in avoiding the above-men- 
tioned ambiguity by an ingenious device to which he was led 
through considerations of S, Eilenberg and S. MacLane’s (1945) 
theory of categories and natural equivalences. 

Whether we interpret the vertices of the graph as ‘‘biological 
functions’? or as the anatomical carriers of those functions, we 
may consider them as components of the biological system. Rosen 
points out that each such component receives either from the en- 
vironment or from other components certain inputs, and emits ei- 
ther to the environment or to other components, certain outputs. 
Thus the property of ‘‘sensitivity to light’? has as input the light 
waves and as output a set of neural discharges. Or we may say 
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just as well that the organ of sensitivity to light has the above- 
mentioned input and output. Thus now the arrows in the graph 
represent not merely relations but inputs and outputs. Those in- 
puts and outputs may be either of a chemical nature, like various 
metabolites, or they may be complex physical or physicochemical 
factors, like light waves, action currents, etc. 

Through the use of the theory of categories, Rosen was led to 
the important concept that every ‘‘biological function’’ or its 
anatomical vehicle represents mapping of certain sets of inputs 
onto certain sets of outputs. Thus in vision a set of wave lengths 
received as input maps on a particular excitation pattern in the 
optic nerve. In digestion, a set of ‘‘admissible’’ proteins maps 
on a set of aminoacids which are the result of digestion. 

Rosen has drawn a number of interesting conclusions which may 
be of biological importance. In all his work hitherto he confines 
himself to one organism. In fact, though not stated too explicitly, 
Rosen practically confines the type of graphs which he has stud- 
ied to unicellular organisms. 

On the other hand, in our studies we emphasized the relational 
similarity between different organisms. We gave expression to the 
fact that any higher organism can be mapped in a many-to-one way 
onto a lower organism with continuous mappings of graphs or of 
topological spaces, or relation-preserving mappings of sets, Inas- 
much as the topological approach is now found to be only one of 
several more general possibilities, we propose here that instead 
of speaking of the principle of biotopological mapping we speak 
of the principle of relational epimorphism in biology. 

Whatever its final precise formulation will be, this principle 
represents one of the most general biological facts and expresses 
the well known unity of the biological world as a whole. It is 
therefore indicated to investigate the relation between Rosen’s 
ideas and the principle of relational epimorphism. That is the 
purpose of this paper. We shall assume the reader’s familiarity 
with all the above-mentioned papers by Rosen and by the author. 

We shall first generalize somewhat Rosen’s theory of the (2, )- 
systems. Before we do that we shall discuss what appears to be 
a rather severe biological restriction used by Rosen. 

Rosen makes the assumption of non-contractibility of inputs. In 
biological terms, this means that in order to function properly 
every component requires ai] the inputs which it receives. The 
failure of a single input, out of many, makes the component in- 
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operative and incapable of emitting its outputs. The assumption 
of non-contractibility is basic for Rosen’s theory. In particular 
the elegant theorem 3 (RI) stands and falls with it. At the same 
time this assumption seems biologically extremely restrictive. 
One has only to think of the common final path in the nervous sys- 
tem. A muscle may be innervated either by pyramidal or extra- 
pyramidal cells. 

It may be argued that though in the case of a double innervation, 
we merely have the same input coming from two different compo- 
nents, it is, however, possible that we have here two kinds of 
inputs which may be different in the patterns of the nervous dis- 
charge. Yet the muscle can function with only one of those in- 
puts. We find similar situations in the humoral system. Of course 
if we restrict ourselves, as Rosen does, to unicellular organisms, 
the above examples are no objections. Yet, in view of the prin- 
ciple of relational epimorphism it would be, so to say, unlikely, if 
not outright impossible, that the hypothesis of non-contractibility 
would hold for unicellular organisms and not hold for multicellular 
ones, 

But do the above examples really contradict the assumption of 
non-contractibility. They do so only if we give the words ‘‘normal 
operation of a component’’ a rather restricted meaning. If we con- 
sider the instantaneous contraction of a muscle, hic et nunc, as 
its normal operation, then the assumption of non-contractibility 
does not hold in this case. But if by normal operation of a muscle 
we mean its over-all normal contribution to the working of the or- 
ganism taken over the whole life span of the latter, then the above 
example in no way contradicts the assumption of non-contracti- 
bility. A harmonious working of an organism may require both 
pyramidal and extrapyramidal innervation. The individual con- 
tractions of a skeletal muscle may be practically identical whether 
innervated pyramidally or extrapyramidally. But the lack of either 
one or the other will not permit the muscle to always perform its 
functions satisfactorily. 

Thus, though the assumption of non-contractibility of inputs may 
lead to some difficulties, we shall preserve it as a useful hypothe- 
sis which may even eventually be elevated to the rank of a general 
principle. 

A problem which looks very much like the one just discussed 
but which is actually of a very different nature concerns the pos- 
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sibility, demonstrated by Rosen (RII, RIII), of bringing the graph 
of an organism into a canonical form in which each component 
emits only one output. 

From a formal mathematical point of view, the procedure is 
Straightforward. Yet the question arises whether it always has a 
biological meaning. A component which produces, say, two differ- 
ent outputs is likely to have in its internal fine structure at least 
two different mechanisms, each responsible for one of the outputs. 
Even if those two different mechanisms partly overlap, we can 
consider them meaningfully as components. Consider, however, 
the following example. A component receives as input a swarm of 
molecules AB. The ‘‘biological function’’ of the component con- 
sists of an enzymatic splitting of AB into A and B, which are 
emitted as two different outputs. What in this case does the for- 
mal reduction to two single output components mean physiologi- 
cally? No matter what the physicochemical mechanism of the 
splitting of AB into A +B may be, A and B are emitted as two out- 
puts, each of which may even go as input to a different component. 

It would seem that such situations will occur relatively less 
frequently than those where a separation into single output compo- 
nents is physiologically meaningful. It seems, however, that the 
reduction to the canonical form may have less value than it looks 
now. Since we are dealing not with mathematical abstractions but 
with biological realities, any mathematical formalism must be con- 
sistently meaningful biologically. This does not seem to be the 
case here. The difficulty, however, does not look insurmountable. 
The introduction of the canonical form removes certain ambiguities 
in the ordinary graph, or block diagram. In considering, however, 
Rosen’s abstract block diagram, it does not seem to offer any 
absolutely essential advantages. 

Until the above-mentioned difficulty is settled, we shall confine 
ourselves to systems in which the canonical form is physiologi- 
cally meaningful for every component. 

Now we come to another restriction in Rosen’s theory of the 
(M,®)-systems which, as we shall see, can be readily removed, 
thus providing a generalization of Rosen’s theory. 

In Rosen’s theory (RI) every component #,, which serves as 
repair mechanism to the associated component M , receives its in- 
puts from some environmental outputs of the (M,®)-system. Con- 
sidering the circumstance that the environmental outputs of any 
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biological system are either mostly waste products of metabolism 
or representative of such phenomena as locomotion, biolumines- 
cence, secretion of protective chemicals, etc., it seems unlikely 
that only environmental outputs would contribute to the reestab- 
lishment of a lost or damaged component. Much more likely, the in- 
puts into the different R,’s may come from outputs of any non- 
terminal component. We shall now show that making this more 
general and biologically much more satisfactory assumption, we 
still preserve all the essential features of Rosen’s theory of 
(N,R)-systems. 

In addition to origins and termini (terminal components) (RI, 
p. 248), we introduce the notion of restituent component, or simply 
restituent. A restituent is a component of the Ji-system such that 
some of its outputs, which in general go to other components, 
serve as inputs to some %, components. Such outputs we shall 
call restituent outputs. Just as in Rosen’s theory a component 
may be both an origin and a terminus (e.g. M, in Figure 1 of RI), 
SO a component may be simultaneously an origin, a terminus, and 
a restituent. Thus the set 0 of all origins in a given block dia- 
gram, the set % of all termini, and the set Q of all restituents have 
in general a non-empty intersection. 

We denote by 6* the set of all restituent outputs. We denote by 
0* the subset of 6*, which serves as input to 2; To each compo- 
nent M, there corresponds the set of outputs ©.* Cc @* which is 
terminated if M; is destroyed or inhibited. We shall call ©* the 
generalized dependent set of M,. 

Now we define a generalized subsystem M* of the ‘l-system 
which is part of the (l,®)-system. By generalized subsystem l* 
we understand a subset of the set of components of St, connected 
in Yl, such that: 

1. The set of restituent outputs of * is a subset of the re- 
stituent outputs of ‘I. 

2. The subset ‘l* receives no inputs from any component of I 
which is not in Jl*. 

Instead of Rosen’s theorem 1 (RI) we now have: 

Theorem 1. If a component M, of an ({l,8)-system is inhibited, 
the possible consequences are: ° 

1. The system fails completely, 

2. There exists a unique entire generalized subsystem Jl* of 
such that §* (together with the corresponding systems, ?., M,€ 
N*) is not affected by the inhibition of M,.. : Z 
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The proof is almost exactly the same as in Rosen’s case, ex- 
cept that for the word ‘‘terminus’’ we substitute everywhere ‘‘re- 
stituent,’’ and for 0; we substitute @*. This is owing to the fact 
that the ‘‘terminality’’ of the components in Rosen’s theory does 
not play any role as such. The importance of terminal components 
lies in the fact that they are also, in Rosen’s theory, restituent 
components, 

In a similar manner we prove Rosen’s theorems 2 and 3. 

The generalized subsystem Yt* which remains functional after 
the inhibition of an M;-shall be called again the related general- 
ized subsystem of M;. The set G* of restitutional outputs of 
tt which are not produced by %,* is called the generalized aug- 
mented set of M.. Just as in Rosen’s case, we see that the neces- 
sary and sufficient condition for a component M, to be re-estab- 
lishable is 6,;* 4 G* = @ or, which is the same, 6; . eo. Gs. 

We now shall investigate the relation Bersean our generalized 
(M,®)-system and the principle of relational epimorphism. Let an 
(M,R)-system represent the primordial organism. We shall con- 
struct from it a system of a higher organism, using the most gen- 
eral method described previously (Rashevsky, 1959a, hereinafter 
referred to as I). . 

We consider the graph of the primordial as being in the canoni- 
cal form, each component producing only one output. 

In this construction we considered subproperties of a given bio- 
logical property (Rashevsky, 1958). If the components of our 
primordial (i,#)-system stand for the various properties, then 
there is no difficulty in speaking of subcomponents, If the com- 
ponents are interpreted as structural vehicles of some properties, 
we again can speak of subcomponents, meaning by ‘‘subcompo- 
nents’’ the structural vehicles of the subproperties. 

The statement ‘‘Component M, emits an output which serves as 
input to component M,’’ describes a binary relation, MRM,, be- 
tween M, and M,. This relation happens to be neither reflexive, 
nor symmetric, nor transitive. From the point of view developed in 
the previous paper (I), the nature of the relation is immaterial for 
our purposes. 

Following now the method of constructing a set which describes 
a higher organism, as discussed in the paper just mentioned, we 
proceed, applying rule (A,) of I, as follows: 


(A) For each component M, of the primordial (i, })-system we 
choose a number n, of subcomponents M,,. If in the primordial 
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(Tl, R)-system, M, emits an output which serves as an input into 
M,, then in the system O(Il,X) which corresponds to the higher 
organism, at least one subcomponent M,, emits an output which 
serves as an input to at least one subcomponent My, g- Further- 
more, for at least one z, an M,,, must emit an output which serves 
as inputs to at least two different M, 2 and M,,,,(B #y). 

As said above, we provisionally adhere to the assumption that 
it is biologically meaningful to represent a system in canonical, 
single output form. We therefore consider that each subcomponent 
M,., emits only one output which is directed to different other sub- 
components ¥, ,.- 

Inasmuch as’ the relation ‘‘M, emits an output which serves as 
input to M,’’ is represented in the graph of the (il, )-system by an 
arrow from M, to M,, we may formulate the above rule A in a purely 
geometric way, namely: 

(A,). If in the primordial (l,%)-system M; is joined to M, by an 
arrow M, — M,, then in the system 6(I,R) at least one subcom- 
ponent M,, is joined by an arrow to at least one subcomponent 
M,,g, SO that M;,, —M,,- Moreover for at least one z a Compo- 
nent M,, must emit arrows to at least two subcomponents My and 
My (6B #y). : 

The #&, components are treated in exactly the same manner. 
Thus if in the primordial (il,R)-system, &, receives an input from 
M 
M 


,» then at least one F,, will receive an input from at least one 


kB° 
We preserve in G(il,) the assumption of non-contractibility of 
all inputs. This has an immediate biological implication. If a 
certain structure in a lower organism is differentiated into several 
substructures, this differentiation does not in general consist of 
a mere addition of new mechanisms in the structure. Each of the 
substructures becomes more complicated, because the number of 
inputs into a component M,, in O(Il,R) is in general greater than 
the number of inputs into ,, and they are all necessary. 

It is readily seen that 6(M,R) is itself an (M,R)-system. We now 
prove the following theorems. 

Theorem 2. Let S(I;*) be the system constructed according to 
rule (A) or (A,) from the generalized related subsystem M,* of M; 
in the primordial. Then 6(%;*) is the generalized related subsys- 
tem in ©(%,R) for the whole set of subcomponents M,,,. 

Proof. First we remark that M;, C G(M) — G(M,;*). If it were con- 
tained in 6(M,*) then, by our construction, M; would be contained 
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in t;*, and thus %,* would not be the generalized related subsys-: 
tem of M, in the primordial. Now suppose the theorem false for 
some M,,,. This means one of two things; either 

(a) there exists at least one M, , in 6(2,;*) which is inhibited by 
the inhibition of any one of the M,,,, or 

(b) there exists at least one M),, in O(N) — G(M,*) which is not 
inhibited by the inhibition of any M,,. 

Consider case (a): An inhibition of the component M, in O(2l;*) 
by the inhibition of M;, can occur only in two ways. One is the 
existence of a directed path from M;, to M;,. This is impossible 
because in that case there would be a directed path from M; to M, 
in the primordial, with M, being in M,*. Then %,* would not be the 
related subsystem of M.. 

The other way is for an output of the generalized dependent set 
6*,, in 6(M) — O(N;*) of some component M,, in O(N) - 6(2,*) 
which is inhibited by inhibition of M;,, to be the input of F, ¢ 
But if this were the case, then in the primordial (J, f)-system, 
there would be a component M, in Sl - 2,* which would be inhibited 
by inhibition of M,; and which would contain in its generalized de- 
pendent set @,;* in J - 2,* an output that is an input into &,. In 
that case, however, ‘;* would not be the related generalized sub- 
system of M,. 

Now consider case (b): It implies that M,,, cannot be reached by 
a directed path from any M;,, nor is there a directed path from 
M;,, to R,,. But this means that there is no directed path between 
M,; and M, in the primordial, nor is #, reached by a directed path 
from M, in the primordial. This again would mean that {],* is not a 
related generalized subsystem of M., 

Thus the inhibition of all M;, leaves G(%l,*) functioning. How- 
ever, an inhibition of a simple M,, may leave more components of 
6(M,#) functional than are contained in G({l,;*). This is so because 
it is possible that there is no directed path from a given M,, to 
M,.,, although there is such a path from M, to M, in the primordial. 
Only if we consider 6(SI,#) to be a maximal set (Rashevsky, 1958) 
will this be possible. Thus we have: 

Theorem 3. If 6(M,#) is maximal, then 6(Mt,*) is the generalized 
related subsystem for every subcomponent M,,. 

Theorem 4. If in the primordial (%l,#)-system M, is re-establish- 
able, then in O(I,R) every M,,, is re-establishable. 

Proof. In the primordial the condition of being re-establishable, 
OF 6 = § states that all inputs into &, come from restituents in 
M,*. By construction, all inputs of any #;, in 6(M,R) come from 
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the restituents of 6(M,*). Hence the condition of re-establishment 
is satisfied. 

Theorem 5. If in the primordial, M, is not re-establishable, then 
none of the M,,, are re-establishable in 6(N,*). 

Proof. By assumption Ff, receives some inputs from It - It,*. 
Then by construction, every f,, receives an input from at least 
one restituent in G(2) - O(Ml,*). Hence, denoting by ©(0, *) the 
set of all inputs received by F;,, in 6(Nl,#) and by 6(G,*) the set ae 
all restitutional outputs of G(Ml)- G(ii,*) we see that G0, a 
‘ay G(G, *) £6. Hence M,, is not re- spetabfiananles 

It may be of interest to generalize the construction in such a 
manner that the number 7,’ of subcomponents #;,, is not the same 
as the number n, of subcomponents M;,. In that case, if n,<7,, 
some F,,,’s will replicate more than one M,,,. If n,’>n, then more 
than one #,, may be needed to replicate one M;,. Biologically 
this is not an implausible situation. 

Thus far we have not made any use of Rosen’s interpretation 
of his (2,R)-systems in terms of categories. Since the theory of 
categories appears to be very promising in relational biology, let 
us investigate briefly the structure of the category, which is im- 
plied by the principle of relational epimorphism. 

The category which represents a higher organism maps in a 
many-to-one way on the category of the primordial. The objects of 
the category are the inputs and outputs. The mappings are the 
biological components. 

The set of objects in the category of a higher organism is thus 
more inclusive than the set of objects in the category of a lower 
one. This reflects the fact that a higher organism is more adapt- 
able, in the sense of being able to utilize more environmental in- 
puts (Rashevsky, 1959b). 

We shall denote by C, the category which represents a higher 
organism, and by C, the category which represents the primordial. 
Thus the functor % ats assigns to each object in C, an object 
in C_ also in a way compares the admissible ayeinnieae of the 
two organisms, as well as the physical and chemical processes 
admissible within the organisms. 

Consider, for definiteness, an origin of the primordial and its 
subcomponents in a higher organism. Each subcomponent maps 
different sets A,,A,,..., of environmental inputs on some of its 
outputs. Thus if the component of the primordial is the general 
sensitivity, then one subcomponent has as domain of its mappings 
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the set of light waves, another—the set of sound vibrations, etc. 
In 6(M,R) the components are interconnected in such a manner 
that many of them (in a maximal ©(%,8) even all of them) even- 
tually affect the same component to which they are connected by a 
directed path. Suppose the first origin subcomponent performs the 
mapping f, = A, —> B,; the second—the mapping jee A, > 85, 
and so forth. Let the component which is reached from all the 
origins map its input onto the set Z of outputs. Since we deal 
here with a category, the law of composition of mappings holds, so 
that if fe©H(A, A’) and ge¢H(A’,A’’), then there is a mapping h = 
gf, such that he H(A, A’’). Hence the set of environmental inputs 
A, is eventually mapped, through a composite mapping fx onto Z; 
thus 7}: A, —>» Z. Similarly for every other A, we have a com- 
posite mapping f;:A,—+ Z. But then it follows that there are 
mappings A, —+ A, for every z and &. In other words the sets of 
mappings H(A,, A,) #0. Now if we take sericusly the whole idea 
of representation of biological functions as mappings, we come to 
the unavoidable conclusion that certain types of sensory inputs 
A, (say, optical) must produce some sensory outputs A, (say, 
acoustical). We have here a general phenomenon of synesthesia 
which follows with necessity from the principle of relational epi- 
morphism coupled with what may be called the principle of repre- 
sentation of organisms by categories. 

Moreover the mappings A, —> A, are not arbitrary but are de- 
termined by the mappings 7. Hence the choice of the sets of dif- 
ferent environmental stimuli of different modalities is not arbitrary. 
There must exist between these sets if not an isomorphisn, at least 
some epimorphism. The structures of the environmental factors 
which an organism can perceive are, in some way, related. Since 
in optical and acoustical factors we have gestalt structures, in- 
variant towards certain ‘‘transpositions,’’ we must expect the pos- 
sibility of olfactory, gustatory, thermal, and touch “‘gestalts.’’ If 
the mappings A, —> A, are 1-1 isomorphisms, then the conclu- 
sion about the relational similarities of processes of vision, hear- 
ing, and other sensations become quite far reaching. 


The author is indebted to Dr. Robert Rosen for a discussion of 


this paper. 
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A DERIVATION OF A ROTE LEARNING CURVE FROM 
THE TOTAL UNCERTAINTY OF A TASK 
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UNIVERSITY OF MICHIGAN 


The derivation of H. D. Landahl’s learning curve (1941, Bull. Math. 
Biophysics, 3, 71-77) from a single information-theoretical assumption 
obtained previously (Rapoport, 1956, Bull. Math. Biophysics, 18, 317-21) 
is extended to obtain the entire family of such curves with the number of 
stimuli M (to each of which one of WN responses is to be associated) as a 
parameter. No additional assumptions are required. The entire family 
thus appears as a function of a single free parameter, k, all other param- 
eters being experimentally determined. The theory is compared with a 
set of experiments involving the learning of artificial languages. An 
alternative quasi-neurological model leading to the same equation is 
offered. 


In a previous paper (Rapoport, 1956), we showed how a learning 
curve, derived by H. D. Landahl (1941) from a neurological model, 
could be derived as a consequence of a single assumption, namely 
that the rate of decrease of uncertainty of response per error is 
constant. That is, we derived an equation formally identical to 
Landahl’s in which our single parameter, &, corresponded to his 
(essentially) single parameter, kd. 

We were not able, however, to show how this parameter de- 
pended on M, the number of independent stimuli with each of which 
one of N responses was to be associated. Hence to get a family 
of curves, formally identical to Landahl’s, we were obliged to 
make the same additional (ad hoc) assumption, namely, k = ne cM 
where 7 and € were now two free parameters. 

Here we will drop this extra assumption and will re-derive our 
learning curve from a slightly different hypothesis. The form of 
each single curve of the family will then remain the same. But the 
family will not contain the experimentally controlled parameter M 
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and the resulting dependence of our parameter (now &’) upon M will 
be similar (though not identical) to Landahl’s. We will then fit 
our family of curves to some experimental data. 

We will assume this time that the rate of increase of the total 
information gained by the learner per error is constant. 

Let A(t) be the total information gained after ¢ stimulus presenta- 
tions, the stimuli being chosen at random from a set of M. 

Then, by our hypothesis 


dh/dw =k, a constant, (1) 


where w is the cumulated number of wrong responses. As in the 
previous paper (Rapoport, 1956), we define the uncertainty of the 
response to a given stimulus in ‘‘nits’’ as U = —log, (1 - dw/dt), 
that is, the negative logarithm of the probability of correct re- 
sponse to any one stimulus.* 

Then the total uncertainty remaining will be U =~—M log, (1- 
dw/dt). If H is the initial total uncertainty of the task, the amount 
of information gained will be 


h=H-U=H+M log, (1 - dw/dt). (2) 
Integrating (1), we get 
h=kw, (3) 


where the constant of integration must vanish, since A(0) = u(0) = 0. 
Eence, substituting (8) into (2), 


kw =H +M log, (1 - dw/dt), (4) 

Exp {(kw — H)/M} = 1- dw/dt, (5) 
du 

dt. (6) 


1— Exp{(kw — HM} 


Integrating (6), we obtain, again considering the initial condition 
u(0) = 0, 


Exp {H/M} 


M 
by ta Sear Sea. REEL ile” 
ee Fe TExptH/M} — 1]-Exp{-kt/M} 4 1— 


“Assuming equi-probability of possible responses, this definition 
coincides with Shannon’s uncertainty (Shannon and Weaver, 1949) A 
‘‘nit’’? is a ‘‘natural bit,” i.e., log, 2 = .693.... 
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If the responses to the M stimuli are independently assigned, 
then the total initial uncertainty of the task must be 


H =M log N. (8) 
Substituting (8) into (7), we obtain 
M N 


w=—lo 


hk °* (N—1)Expi-ki/M}4 1° 


(9) 


which is Landahl’s equation (11), where his parameter kd coincides 
with our parameter &’= k&/M (Landahl, 1941). In other words in- 
stead of Landahl’s assumed exponential dependence upon M (kb = 
ne S™), we obtain a hyperbolic relation as a consequence of our 
hypothesis. When only a few curves are available (in Landahl’s 
experiment—three), the distinction between the two relations can- 
not be made with confidence. 

If the responses are not independently assigned, we cannot 
assume (8). However, the quantity e”/™ is still a constant in a 
given experiment. If we designate this constant by N’, our equa- 
tion (7) is still formally identical with Landahl’s equation (11). 
However N’ now becomes the ‘‘effective’’ N, which will, in general, 
be smaller than the actual number of possible responses to each 
stimulus. The decrease is a result of the redundancy introduced 
by non-independent assignments of responses to stimuli. 


Application to an Experiment. Several groups of subjects (10- 
50 in each group) were taught several artificial languages re- 
spectively by associating irregular geometric figures with their 
names in the languages. The referents were the same in all cases, 
namely, a set of figures of five different shapes and five different 
colors. Out of a set of 25 possible combinations, 20 were used as 
the referents, the remaining five as test stimuli wherever generali- 
zation of the naming principle could be made. In the course of the 
training, each of the 20 referents chosen in random order was 
presented to the subject, who was required to make a verbal 
response—a nonsense syllable or a pair of nonsense syllables 
depending on what language was being learned. Only a completely 
correct response (i.e., correct pronunciation) was counted as 
correct. All others were counted as errors. Whether the subject’s 
response was correct or not, the experimenter countered with the 
correct response. Each session lasted about an hour. Successive 
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sessions were held on consecutive days. It took about 4 to 30 
sessions to complete the learning to criterion (three perfect con- 
secutive scores on the 20 responses in random order), depending 
on the language and, of course, on the subject. The languages 
were structured as follows: 

L,: With each of the 20 figures a single syllable was associated. 

L,: With each of the 20 figures a pair of distinct syllables were 
associated, no single syllable being associated with more than 
one figure. 

L,: With each of the five shapes and with each of the five colors 
a single syllable was associated respectively. The syllables 
were combined in pairs to name each figure by shape and color. 
The shape syllable always came first. 

L,: The vocabulary of this language was the same as that of L, 
(10 syllables), but there was no systematic correspondence be- 
tween the single syllables and the aspects (color and shape) of 
the referents. Thus to each referent there corresponded a unique 
ordered pair of syllables, but there were inconsistencies in corre- 
spondences between single syllables and the aspects of the 
reference figures. 

To test our theory, we turn our attention first to the asymptotic 
values of w. Setting ¢ = in equation (7) we get the theoretical 
asymptotic value, 


w(o) =W = H/k. (10) 


If now & is indeed independent of the complexity of the task, then 
H/W ought to be constant for all languages or, at least, its fluc- 
tuations ought to be attributable to population sample fluctuations. 

To test this hypothesis, we need to calculate the total un- 
certainty H in each language. Note that the assignments of re- 
sponses to stimuli are not made independently, except in L,. In 
L,, for example, having assigned a syllable to one figure, we have 
only 19 syllables to choose from to assign one to the next figure. 
Obviously the total number of assignments is 20! Therefore the 
total uncertainty in L, is 


H(L,) = log, 20! = 42.3. (11) 


The task of learning L, was given to 20 subjects. The average 
total number of errors per subject was observed to be 304. Hence, 
the value of & was calculated as 42.3/304 = 0.139. 

Turning to L, and taking the syllable as our ‘‘unit,’? we find 
that the amount of uncertainty in that task is just twice that of Ly 


ROTE LEARNING CURVE 89 


or 84.6 nits. We should add another bit (=0.7 nits) to represent 
the uncertainty of which group of 20 syllables represents the 
repertoire from which the first syllable of the two is chosen, 
since the two syllables must be pronounced in the correct order. 
Taking 85.3 nits as our total uncertainty and observing the total 
average (10 subjects) number of errors to be W = 540, we obtain 
k(L,) = 85.3/540 = 0.154, which agrees fairly well with the values 
previously obtained. 

In L, we had data on 50 subjects. The average W was observed 
to be 85.8. To calculate H, note that the units are now syllables 
to be associated with the ‘‘aspects’’ of the referents, namely the 
five shapes and five colors. Therefore, the amount of uncertainty 
is now 


H(L,) = 2 log, 5! + 0.7 = 10.3. (12) 


Hence, &(L,) = 10.3/85.8 = 0.120. We see that this value fluc- 
tuates, but so far not excessively. Anticipating statistical analy- 
sis, we conjecture that the fluctuations of & observed do not 
exceed those expected in samples of populations of this size. 
Finally, turning to L,, we note that here the responses are as- 
signed independently, each of the two syllables being selected at 
random (with replacement) from its own repertoire of 5. Therefore 


H(L,) =2 log, 5 x 20+ 0.7 = 65.1. (13) 


Observing W = 720 (10 subjects), we obtain &(L,) = 65.1/720 = 
0.09. This is the largest deviation of & we have observed. It can 
perhaps be attributed to the confusion in the language due to the 
non-unique associations between the syllables and the shape-color 
aspects of the referents. 

We now wish to compare the theoretically predicted learning 
curves with our data. We can proceed in two ways. 

1. Determine & from the data on one of the languages. Having 
thus fixed the %, use it in all the experiments to calculate the 
value of H(H = kW), which will then be the free parameter in each 
curve. 

2. Calculate the value of H in advance, as we have done above, 
and use it and W to determine &, which will then be the free 
parameter. 

Comparisons of theory and experiment are shown in Figures 1-4. 
Curves I result if & is kept fixed throughout at 0.139 obtained 
from data on L,, and H is treated as a free parameter. fe Hix, lig, 
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FicurE 1. Comparison of theory and experiment for Ly. Number of 
trials (horizontal) vs. cumulated errors curve calculated from equation 
(7) fitted with k as a free parameter against observed values (circles). 
Here curves IJ and II are identical (see text). 


and L,, M = 20 (twenty referents), but in L,, M = 10 (the referents 
are five shapes and five colors). Curves II result if H is treated as 
an experimentally controlled parameter, calculated for each lan- 
guage as shown above, and £ is treated as a free parameter. In 
L,, L,, and L,, the two curves are barely distinguishable; but in 
L, the discrepancy is considerable. Since in the latter case curve 
II gives definitely the better fit, we prefer to use method 2, which 
incidentally is also more attractive from the point of view of 
applying information theory to the learning model, since in this 
method the ‘‘information’’ of a language is theoretically calculated. 

We note, further, that the theory fits the curves obtained in L, 
and L, very well with relatively slight fluctuation of the parameter. 
It also fits L, very well but with a considerable fluctuation of 
the parameter & in the direction expected because of the ‘‘inter- 
ference’’ inherent in that language. 

In the case of L., however, there seems to be a discrepancy 
between the theoretical curve and the data. 

To discuss the possible reasons for the discrepancy, it will be 
useful to prove the following: 


te) re) 
Lemma. (=) m0). (=) 2 0. 
OH) » ak/ y 
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Since H = kW, it is sufficient to prove the first of these in- 
equalities. 
Proof. Substituting & = H/W into (7), we obtain 


MW 
te Ws =yynlog (ene Meret) eg te 4 4], (14) 


Differentiating with respect to ¢, we get 


ae (eH /M _ 4) e7Ht/WM 


Ae Ey Se er a 15 
dt (eH/M_4)e"Ht/WM 4° oe 


For ¢=0, (15) becomes 1- Exp {-H/M}. Let now w, = w(H,), 
Ww, =W(H,), where H, > H,. Then 


dw des) eS = 

eee he pete tat Ms 5 tal M 9) 16 

e at : : a 
t=0 

600 te 


500 


400 


pag pug 
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FicurE 2. Comparison of theory and experiment for Lg. Curve I is 
obtained by taking the value of & from L, and treating H as a free param- 
eter; curve II by calculating H theoretically and treating k as a free 
parameter. The discrepancy between theoretical curve and data is 
discussed in the text. The discrepancy between the values of W is due 
to rounding-off error. 
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FicureE 3. Comparison of theory and experiment for L3. Here M = 10. 


In other words, w, increases faster than w, at ¢=0. Further, if 
W is fixed, w,(0) = w,(~)=W, and of course w,(0) = w,(0) = 0. 
Our lemma will, therefore, be proved if we can show that w, - w, 
has a single maximum, i.e., (dw,/dt —- dw,/dt) has a single posi- 
tive root. Setting this.expression equal to zero, we obtain from 
(15), after simplifications 


(elale 211) ea t/WM _ (erum e 1)e-# it/WM _Q, (17) 
La 
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Figure 4. Comparison of theory and experiment for L4. Here the dif- 
ference between curves I and II is considerable. Preference, therefore, 
is given to calculating 4 theoretically and taking a lower value of k. 
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But this is of the form 


ae *'—~ be ¥' = ae ** 3 ~ 5 tem =0, (18) 
a 
where a> 6, x>y. The expression in the brackets vanishes for 
exactly one positive value of ¢, and its coefficient never; so our 
lemma is proved. 

The result can now be applied to our theory as follows. If w(0) 
is pegged at 0, and w(~) at some observed asymptotic value W, 
then if H (and consequently &) increases, the entire curve shifts 
upward and vice versa. 

Cbserve now the discrepancy between theoretical and experi- 
mental curves in L,. This discrepancy can be reduced (perhaps 
practically eliminated) if the value of & (and of H with it) is taken 
lower than 0.154. We note that this value is the highest value of 
k obtained. Consequently a lower value will be more in accord 
with the other values of & obtained in the other experiments. But 
is there a justification for reducing H, which, after all, was calcu- 
lated by the same method as in the other experiments, where the 
agreement was good? It can be argued that specifically in the 
case of L,, H can indeed be taken lower than we have taken it. 
The assumption in calculating H(L,) was that the stngle syllables 
were the units of which the repertoire is composed. But in this 
language a pair of such syllables is inseparable. Hence, logically 
speaking, a pair of syllables could be taken as a unit. If this 
were done, L, would have exactly the same H as L,, namely 42.3 
nits. In that case, it can be shown that the theoretical curve for 
L, will lie entirely below the data. The conjecture is that the 
pairs of syllables ‘‘coalesce’’ but not completely. An intermediate 
value of H between 42.3 and 85.3 should give a fit of accuracy 
comparable to that obtained in the other languages. 

Note that the argument about the inseparable syllables does not 
apply either to L, or to Ly, the other two-syllable languages. 


Theoretical Relation to Landahl’s model. It might be construed 
that our basic assumption, dh/dw =k implies that the subject 
‘fearns’’ (i.e., gains information) only from his wrong responses, 
in contradiction to the assumption underlying Landahl’s equation 
(11), which states that increments of excitation at the synapse 
leading to the correct response result only from the subject’s 
correct responses. 
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In view of the formal identity of the equations derived from each 
of these assumptions, this seems anomalous. However, it is not 
true that our assumption implies that the subject gains information 
only from wrong responses. To show this let us calculate the 
rate of change of information with respect to correct responses c. 
This is 

dh/de = (dh/dw) (dw/dc) = kdw/de. (19 ) 


Now 


»/dt 
enlelalere (20) 
1 - dw/dt 
Combining (15), (10), (20), and (19) we get 
es Kent aye eX, (21) 


dc 


In other words, our model implies that the rate of gain of in- 
formation per correct response is an exponentially decaying func- 
tion of ¢, the number of trials. This result is intuitively plausible. 
At the start of the learning process, when very few correct re- 
sponses have become fixed, the rate of gain of information per 
correct response is large because of the possibility that some 
responses, correct for the first time, will become fixed. But at 
the end of the process where most of the correct responses are 
already fixed, the gain of information from correct responses is 
small. In particular no information is gained from a correct re- 
sponse which had already been fixed. 

Another way of looking at the model is to find the expression 
for dh/dt. We note that dh/dw = (dh/dt)/(dw/dt) = (dh/dt)/(1 - P_), 
where P. is the probability of correct response. If dh/dw =k, a 
constant, we have dh/dt = k(1— P.), i.e., the information increases 
proportionately to the number of trials and the amount ‘‘left to 
learn.’’ This seems reasonable for the case of 100% prompting, 
which is the case in our experiment. 


Relation to Estes-Bush-Mosteller models. The simplest form of 
the type of rate-learning models proposed by W. K. Fstes (1950), 
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R. R. Bush and F. Mosteller (1955) is expressed by the following 
recursive formula: 


Piyyp =%, + (1-G)A. (22) 


Eere p is the probability of some response, and 2d are con- 
stants. If p represents the probability of the correct response, 
and if the process ends in perfect learning, X must equal unity. 
Hence, in this case, subtracting p, from both sides, 


Peat ~P, = (1 - 4%) -(1- &)p. (23) 


The emphasis in the theoretical development of the authors 
mentioned is on stochastic aspects of the learning process. That 
is, the center of attention is on the frequency distributions of the 
responses on the successive stimulus presentations. 

Cur model, on the other hand, is deterministic. We take for our 
time dependent variable, the average number of responses (say 
over a large population) at each stimulus presentation. 

Passing from the difference to a differential equation, we obtain 
from (23) 


dp 
— =(1- “)(1-p). 24 
mrt )(1- p) (24) 
Integrating, 


pete ts Jay Seammieree (25) 


where p, is the initial probability of the response in question. 
Now in our notation p (the probability of correct response is 
dc/dt, and 1- pis dw/dt. Substituting into (25), we obtain 
we = Ae Sere, (26) 
dt 
where A =dw/dt at ¢ =0 is the initial probability of the wrong 
response. Integrating once again, we obtain the expression for w 
based on the deterministic version of Fstes-Push-Mosteller model, 
namely 


aed (eae eg ae), (27) 
1-a 
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Since A is experimentally controlled, equation (24) contains one 
free parameter 4, as does our equation (7) and Landahl’s equiva- 
lent equation (11). The parameter can be fixed by observing the 
asymptotic value of w in the learning curve. All attempts to fit 
the experimental data described above by equation (27) show a 
large discrepancy, all in the same direction: the theoretical curves 
are all considerably below the experimental. We suspect that our 
model and Landahl’s, which are formally equivalent although 
couched in different terms, are more suitable for the learning 
curves of human subjects, since they take into consideration (as 
our analysis above has shown) the characteristic fixating of the 
responses as they are learned in the process. 

Cf course in the theory of stochastic models, equation (27) 
applies only when there are only two possible responses, and 
moreover the effect upon the probability of the correct response 
is the same whatever response is given. Because of the very 
special assumptions which underlie equation (27), the failure of 
this equation to predict our learning curves in no way invalidates 
the previous models. It only shows the inapplicability of the very 
special assumptions in this case. 

It will nevertheless be instructive to compare the equation de- 
rived from our model, which is analogous to equation (27), that is, 
the expression of dp/dt in terms of p, with (27). 

From (7), we get by differentiation 


dus (eH /M _ 4) e7kt/M 


dt (e#7M — le tM AS (28) 


In the notation of Bush and Mosteller, dw/dt is 1-p. Also 
ef/M is (1- A) '. Hence, translating (28) into that notation, 


pw (Bertha yyc!, (29) 
where B = A/(1~ A). Differentiating, 

dp Bhas em 

=, 30 

dt M(Be*™ , 4)3 em 
Substituting (29) into (30), we obtain 


k 
ae ye - P) (31) 
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But (81) is the well-known logistic equation. It represents the 
simplest case of a contagion process. Based on this analogy, one 
can now give our model a quasi-neurological interpretation, similar 
to the one given by Bush and Mosteller (op. cit., Chap. 2). 

Suppose the probability of correct response is given by the 
fraction of ‘‘properly conditioned’’ neural elements, i.e., a correct 
response occurs if the response happens to be channeled into a 
“properly conditioned’’ pathway, but not otherwise. Then our 
equation (31) says that the properly conditioned pathways behave 
like ‘‘infected’’ individuals who infect others until every pathway 
is ‘properly conditioned.’’ 

Equation (24), on the other hand, derived from the simplest as- 
sumption of the Fstes-Bush-Mosteller type, says that the ‘‘infection”’ 
strikes randomly from the outside, i.e., the fraction of the non- 
infected is exponentially decaying. Our conjecture therefore that 
equation (7), or equivalently, (31), is more characteristic of human 
subjects may well be a conjecture about a higher interdependency 
among conditioned pathways in the human brain. 


This work was supported in part by a research grant from the 
Rackham Fund of the Horace E. Rackham School of Graduate 
Studies, University of Michigan, and in part by a grant-in-aid by 
the International Society for General Semantics. 
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The rapid growth of quantum mechanics in the twentieth century has in 
no small part been due to the advanced state of classical mechanics. 
The introduction of new mathematical concepts to biology and medicine, 
on the other hand, is difficult because of the absence of a strong mathe- 
matical foundation. As a result, in many instances, new theoretical ap- 
proaches to biological material have depended upon the use of mechanical 
and electrical analogues. As this paper shows, from a solution to a 
mechanical problem by means of quantum mechanics comes a potentially 
useful approach to a common medical problem. 


Theoretical Results. In classical mechanics, the motion of a 
particle is strictly predictable from its initial behavior and its en- 
vironment. Thus, a particle of initial energy E will invariably 
pass over a single-step potential barrier of height V, < E. How- 
ever, according to quantum mechanics, reflection of the particle 
may occur even if E >V,. Now, T is defined as the probability 
that a particle in a given state will pass by a given barrier, while 
R =1-T is the probability of reflection by the barrier. It can be 
shown (Rojansky, 1938) that for a particle of energy FE impinging 
on a single-step barrier of height V,, 


4J1-V,/E 


~ a STAVE - 


As the table and figure show, the transmission is imperfect unless 
FE is much greater than V,. 

Stimulation of a biological organism ordinarily leads either to a 
continuous or an all-or-none response. When the change in be- 
havior is a continuous, reversible function of the modifying stimu- 
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lus, the mathematical approach is straightforward. For example, a 
normal distribution is apt to represent closely the change in the 
pulse rate or the weight of a large number of rats produced by a 
given dose of a drug. Furthermore, a logarithmic relation is apt to 
exist between the means of the normal curves and the dose. In 
contrast, the mathematical representation of an all-or-none response 
to stimulation is presently quite inadequate. Typical examples of 
this type of behavior are the appearance of a rash and the death of 
an animal. The all-or-none response is a distinctly different event 
from the gradual changes leading to it and is not readily reversible 
by decreasing the stimulation. Its appearance represents a dis- 
continuity in the mathematical sense. Evaluation of this pattern 
requires administration of a constant dose of the drug to a large 
group of similar animals and observation of the percentage which 
responds in the all-or-none fashion. As the number of animals in- 
creases, this percentage approaches the probability of the response 
for the given dose. Obviously, with only two possible states—the 
initial and final ones—a normal distribution does not exist. 

The relation between the probability of the response and the 
dose is of interest. For small doses, the response does not at all 
occur. Eventually, as the dosage increases, the desired effect be- 
gins to appear in a small percentage of the test animals. A dose 
only slightly greater would be expected to cause the end response 
in most of the animals, since they are so much alike. The proba- 
bility of response would asymptotically approach one as the dosage 
increases further. This pattern resembles the results in Table 1 and 
Figure 1 closely. The resemblance suggests application of the quan- 
tum mechanical model by analogical reasoning. Thus, each individ- 
ual organism is assumed to offer a constant barrier V, per unit mass to 
the appearance of the end response. For unknown reasons, stimuli 
greater than the barrier do not always produce the final state just 
as, with E >V), the particle fails to go over the barrier in every 
case. The determination of the best value of V, would depend upon 
curve-fitting techniques. Ordinarily this could best be done graphi- 
cally by means of a plot of the observed probability as a function 
of the dose. Superimposed upon this curve would be the curve pro- 
duced by equation 1, plotted with various values of E chosen as 
V,- The value of V, would be determined from the theoretical 
curve with the best fit. This value would represent the constant 
intrinsic barrier to attainment of the end-response. 


QUANTUM MECHANICS AND MEDICINE 101 


TABLE 1 


Values of 7 as a function of E/ oe 


ELV z 
1.00 0.00 
1.03 0.50 
soit 0.73 
2.00 0.98 
o0 1.00 
~ 
1009 
fe) 1 2 P=E/V, 


FIGURE 1. A graph of T as a function of E/V,. 


Discussion. In a study of all-or-none behavior, only two dis- 
continuous alternatives exist for each animal—either to respond 
or not to respond. Despite this sharp dichotomy, the percentage of 
animals who respond in all-or-none fashion to a given amount of a 
drug may vary continuously as the dose changes. Unfortunately, a 
paucity of this type of dose-response curve exists in the literature, 
for unknown reasons. Instead, all-or-none behavior is characterized 
by such arbitrary magnitudes as the LD 50—by definition, the 
dose of a drug which causes the death of fifty per cent of the 
animals who receive it. The LD 50 tells nothing about the re- 
sponses to other doses of the drug. 

It is possible that the suggested quantum mechanical approach 
is merely a convention. Physics, as has been shown by Poincaré, 
has many such conventions, which cannot, because of their very 
nature, be either proven or disproven experimentally. A conven- 
tion is used because it simplifies interpretation of related physical 
phenomena. Thus, empirical use of the value V, to characterize a 
species would, at least, indicate that Table 1 represents the shape 
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of the response pattern. Where the LD 50 represents merely the 
response to a single dose, V, would represent a whole pattern of 
response. Projected experimental studies will show whether or 
not Table 1 is actually as closely representative of all-or-none be- 
havior as it theoretically appears to be. The practical usefulness 
of the quantum mechanical approach as a convention will depend 
upon the results of these experimental studies. 

The possibility also exists that the concept of constant species 
resistance to an all-or-none stimulus is intrinsically true. In the 
classical viewpoint, no exact statement can be made about the 
intrinsic resistance of an organism to a given stimulus. Thus, the 
fact that a group of animals have the same all-or-none response to 
the same dose of the drug does not mean that the animals have the 
same intrinsic resistance. Rather, since each animal could not, 
for technical reasons, receive an absolutely identical dose of the 
drug in an exactly identical way, the sameness of the end-response 
could only be due to a cancelling out of variable factors. In the 
classical viewpoint, as the stimulus approaches exact uniformity 
in increasingly uniform-animals, the range of dosage producing a 
percentage response of more than zero but less than 100 per cent 
would become narrower and narrower. As experimental fastidious- 
ness increases in the quantum mechanical approach, the response 
would, instead, tend to be more consistently close to the curve of 
Table 1. Experimental validity of the quantum mechanical concept 
would imply that perfection in biological experimentation is not 
possible and that experimental imperfection rather than variability 
in intrinsic resistance causes variability in all-or-none response. 
This conclusion would be similar to the Heisenberg Uncertainty 
Principle in quantum mechanics, which states that, at atomic 
levels, the process of measurement significantly distorts the quantity 
measured, The possibilities raised suggest an experimental study 
whose results cannot be anticipated at this time. 
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A reappraisal of recently proposed definitions and criteria of active 
transport in terms of experimentally accessible parameters leads to the 
conclusion that it is in principle impossible to give a rigorous quanti- 
tative definition of the active component of the flux of a specific molecu- 
lar species across a membrane without prior knowledge of the mechanism. 
The attempted distinctions between various types of ‘‘non-passive’’ 
transport (coupled, forced, facilitated, etc.) are thus of necessity oper- 
ationally ambiguous. The essential equivalence of formulations based 
on classical thermodynamics and on thermodynamics of irreversible 
processes is pointed out, and the significance of Curie’s theorem to any 
theoretical formulation is discussed. 


Despite numerous experimental and theoretical investigations in 
the field, the concept of active transport has thus far eluded a 
rigorous definition, i.e., one which is both conceptually unam- 
biguous and at least in principle subject to direct experimental 
test. The limitations of T. Rosenberg’s (1948) definition of active 
transport as transport against an electrochemical potential gradient 
and of its modification given by H. H. Ussing (1949) are now well 
recognized (cf. Rosenberg, 1954; Ussing, 1955; Glynn, 1956; Jar- 
detzky, 1957; Patlak, 1956), but no operationally meaningful alterna- 
tive has been put forth expressing the intuitive notion of transport 
by a specific mechanism for which chemical energy is expended. 

Since mathematical treatments of phenomena associated with the 
term ‘‘active’’ transport (frequently containing elaborate schemes 
for the analysis of experimental data; cf. Patlak, 1956; Scheer, 
1958; and others) periodically appear in the literature, it seems 
appropriate to examine the question as to what minimal experi- 
mental information must be available in order to permit a meaningful 
formal analysis of the phenomenon. The measurements most readily 
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available are those of a flux of a given component and, at least in 
some cases, of the associated thermodynamic forces. The question 
then is whether it is possible on the basis of such measurements 
(a) to ascertain the occurrence of ‘‘active’’ transport and (b) to 
calculate the magnitude of the ‘‘active’’ component to the flux. 
The problem is complicated by the fact that in formal treatments 
the term ‘‘active’’ can be and has been used with either of two 
quite distinct connotations, i.e., (a) that the observed flux of a 
species j cannot be fully accounted for in terms of the thermody- 
namic forces known to act in the system, therefore requiring a 
source of energy to account for it or (b) that the observed flux of 
j is fully or in part the primary result of a specific metabolic mech- 
anism, rather than an effect secondary to the transport of some 
other species 7. Thus there are at least four distinguishable ques- 
tions which one may wish to have answered from flux and force 
measurements: (1) Does an energy requirement exist to account for 
a portion of the flux of a species j? (2) What is the magnitude of 
this energy requirement? (3) Is there a specific metabolic mech- 
anism by which the tramsport of j is accomplished? (4) If so, what 
is the magnitude of the flux component to be accounted for by this 
mechanism? 

Much if not most of the existing confusion concerning definition 
of active transport has been caused by attempts to infer from the 
answer to one of the questions an answer to one or all of the 
others. The questions and their answers, however, encompass only 
partially overlapping classes of phenomena (cf. Jardetsky, 1957) 
and therefore cannot be related to one another without additional 
information. 

A theory in which the energy-requirement, rather than additional 
knowledge of the structure of the system, is used exclusively 
as a basis for classifying transport processes in living systems 
has presented an attractive possibility (Rosenberg, 1948, 1954; 
Ussing, 1949, 1955). However, the obstacles to the general and 
unambiguous use of such a theory are fundamental and insurmount- 
able for two reasons: 

(1) A definition stated in terms of the energy requirement derived 
from fluxes and forces without further qualification necessarily ex- 
cludes a large and biologically important group of phenomena in 
which metabolically linked transport occurs in the same direction 
as that of the relevant thermodynamic potential gradient, i.e., the 
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energy may be supplied, even though a requirement for it is not 
obvious. 

(2) Such a definition just as necessarily includes a large class 
of phenomena of only minor relevance to specific transport mech- 
anisms, i.e., those in which the energy required for transport 
against a thermodynamic gradient is supplied in a manner not 
specifically characteristic of a biological system—e.g. by the 
maintenance of compensating electrical, temperature, or pressure 
gradients; by trapping the transported species on one side of the 
membrane, through chemical binding, clathrate formation and the 
like; or by compensatory exchange of one component for another 
(e.g. to preserve electrical neutrality). 

The fundamental difficulties involved are not avoided by a re- 
cent attempt of B. T. Scheer (1958) to give a new classification of 
transport processes on the basis of an apparent modification of the 
definition previously formulated by Rosenberg (1948, 1954). Start- 
ing from the general formula for the rate of entropy production as 
given in the theory of irreversible processes (De Groot, 1952; 
Prigogine, 1955), i.e., 

dAS il 

ee ep nee (1) 

i 

(where S is the entropy, T the absolute temperature, ¢ the time, 
J, the generalized flux of, and X, the generalized force acting on 
the 7*® component), he defines active transport of component j as 
transport under the condition that the product of the flux of com- 
ponent j by the generalized force acting upon it (across the mem- 
brane) is negative, i.e., J,X; “ah i 

It is easy to show that the proposed statement is not truly a 
modification of Rosenberg’s definition, but is in essence equivalent 
to it. Since both J, and X, are vectors in the same space, the 
positive direction is uniquely defined for both by the relation 


J,= Se (2) 
i 


where Li denotes a coupling coefficient (cf. De Groot, 1952; 
Prigogine, 1955). The statement that J, xX, is negative means simply 
that the flux is occurring in the direction opposite to the direction 


of the force. The force in question being the negative of a thermo- 


106 OLEG JARDETZKY 


dynamic potential gradient, the statement reduces to that of Rosen- 
berg, if X is understood to stand for the negative of the electro- 
chemical potential gradient. Similarly, the mathematical criteria 
formulated by Scheer are equivalent to the criterion given previ- 
ously (Jardetzky, 1957). It necessarily follows that the proposed 
definition is subject to all of the limitations and ambiguities which 
have already been discussed (Rosenberg, 1954; Glynn, 1956; Patlak, 
1956; Jardetzky, 1957).* Thus the existing mathematical formula- 
tions of active transport (cf. Rosenberg, 1954; Ussing, 1949; Jar- 
detzky, 1957; Scheer, 1958) provide criteria solely for establishing 
that the flux of a given species j cannot zn toto or in part be ac- 
counted for by the action of such thermodynamic forces as are 
known to act in the system and have been taken into considera- 
tion. It might appear that this would lead to a perfectly rigorous 
quantitative definition of the active component of the flux (and 
hence of the corresponding energy requirement) as the component 
unaccounted for by al/ the known thermodynamic forces. However, 
this would imply that the ‘‘active’’ component of the flux could be 
accounted for by invoking a new type of force, i.e., a force other 
than a thermodynamic potential gradient, which arises as a result 
of metabolic reactions. This possibility can be categorically ex- 
cluded on the basis of Curie’s theorem (Curie, 1908; also cf. De 
Groot, 1952, p. 77).** No chemical reaction or sequence of chemi- 
cal reactions can impart directional properties to the motion of a 
chemical species except by altering the spacial gradient of a 
thermodynamic potential. 

This means that it must be possible either (a) to identify the 
active transport term for a component 7 with a member of the sum in 
equation (2), or (b) account for it by the flux of a distinct com- 
ponent, %, which is derived from j.*** Thus, clearly, a categoriza- 


*The criterion proposed in eq. 3 (Jardetzky, 1957) is of course subject 
to the same ambiguities. Scheer (1958) is however in error, stating that 
it in itself is not wholly reliable; it is, in fact, a rigorous restatement of 
Rosenberg’s definition. 

**The applicability of Curie’s theorem to this problem has been inde- 
pendently pointed out by F. M. Snell (personal communication, to be 
published in ‘Disturbance of Cellular Metabolism in the Premature In- 
fant,’’ Report of the 32nd Ross Conference on Pediatric Research, Colum- 
bus, 1959). 

***For purposes of formal description any alteration of structure or 
configuration associated with a change in energy may be considered as 
the formation of a new component. 
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tion of the flux into ‘‘active’’ and ‘‘passive’’ based solely on the 
possibility of accounting for it by known vectorial thermodynamic 
forces is meaningful only as long as some of the vectorial forces 
remain unknown and may be presumed to arise from metabolic reac- 
tions. It becomes impossible as soon as all the forces become 
known. This, however, does not mean that any distinction what- 
soever between ‘‘active’’ and ‘‘passive’’ components of the flux 
vanishes in a completely described system. It only means that the 
information required to make this distinction possible is not con- 
tained in the transport equations per se. 

The essential distinction between ‘‘active’’ and ‘‘passive’’ 
terms lies in the de facto dependence of the former on a supply of 
metabolic energy. This categorization does not necessarily coin- 
cide with a classification based on differences in the mode of 
transfer. It is for example entirely conceivable that a given spe- 
cies would be transported jointly with a ‘‘carrier,’’ part of the 
energy for this process being derived from metabolic reactions and 
part from other sources. Nor, as already stated, does it coincide 
with a classification based on the apparent existence of an energy 
requirement. It thus hinges on specific prior knowledge of the 
coupling of metabolic reactions to transport processes. The theory 
of such coupling and the resulting detailed formal analysis of 
transport phenomena in living systems will be given in a separate 
communication jointly with F. M. Snell. 

Some important conclusions, the general validity of which does 
not depend on the knowledge of any specific mechanism, may never- 
theless be drawn on the basis of the above analysis: 

Since the information permitting one to categorize transport 
processes into those dependent on and those independent of meta- 
bolic reactions is not contained in the formal description of trans- 
port in terms of fluxes and forces, no criterion based solely on flux 
and force measurements can yield information on the metabolic de- 
pendence of a given transport process. The ambiguity of the infer- 
ences which may be drawn from such classifications as have been 
proposed so far (‘‘coupled,’’ ‘‘forced,’’ ‘‘facilitated,’’ etc.; for de- 
tails cf. Scheer, 1958; Patlak, 1956) seriously restricts their sig- 
nificance. 

With the success of calculating the magnitude of the ‘‘active’’ 
(i.e., metabolically dependent) term depending entirely on prior 
and independent knowledge of the transport mechanisms involved, 
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the question must still be answered as to whether it is possible in 
the absence of such knowledge to establish unambiguously that a 
given species has been actively transported and perhaps to esti- 
mate a limiting value for the extent of such transport. Two kinds 
of assumptions have figured prominently in attempts to arrive at 
such a limiting value: 

(i) That the transport observed after inhibition of metabolism is 
a valid measure of the passive flux or, alternatively, 

(ii) That ‘‘active’’ transport occurs in only one direction, coupled 
with 

(iii) That the ‘‘passive’’ permeability coefficients (i.e., those 
describing ‘‘free diffusion’’) in the two opposite directions are 
equal. 

The impossibility of rigorously distinguishing between an active 
and a passive component of the total flux by using inhibitors of 
metabolism has been pointed out in a previous paper (Jardetzky, 
1957). The assumption of equal ‘‘passive’’ permeability coeffi- 
cients in the two opposite directions under conditions not involv- 
ing thermodynamic equilibrium (which is to say under all conditions 
in which transport occurs), is equally tenuous, quite apart from the 
fact that it neglects all cross-coupling terms of equation (2). Again 
it must be recalled that permeability coefficients (a better designa- 
tion would be apparent permeability coefficients) of a biological 
membrane are not physical constants, and their dependence on the 
thermodynamic potentials of the surrounding media is not predict- 
able a priori (cf. Jardetzky, 1957). From the requirement that any 
net flux be in the direction of the activity gradient, it merely fol- 
lows that the only relation between the coefficients which must be 
satisfied under given non-equilibrium conditions is the following: 


<< 
ij x; 
aE See (3) 
L pe 


=> 
where L,, is the permeability coefficient in the direction of the flux 
and X, the thermodynamic activity on the side from which the flux 
Gite DaloRs In the absence of equilibrium and within this limit, 
LS L,,. 

eas 

This relation holds equally well for non-living as for living sys- 
tems and therefore does not imply any active, i.e., metabolic energy 
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dependent, component of transport. Restricted diffusion (Hodgkin 
and Keynes, 1955) amply illustrates the point. 

The conditions under which it may be possible to conclude that 
a flux component of species j has resulted from metabolic reac- 
tions and to obtain a limiting value for this component can be as- 
sessed on the basis of the following considerations. If the system 
is not completely described, equation (2) takes the form 


=>) L..X,+5;*, (4) 


where the sum is taken only over the m known thermodynamic forces, 
and J * represents the unaccounted-for component of the flux. Two 
assumptions are inherent in the acceptance of J; as a measure of 
the ‘‘active’’ flux, J ,: (1) That none of the known forces X, derive 
from metabolic reactions and (2) that no forces other than those 
resulting from metabolic reactions contribute to J j*. If the second 
of the assumptions holds, but not the first Jx¢<J,. If the first 
does, but not the second J;* >J,. Ifboth do, J = J, and if neither, 
J;* zJ,. Since only J; and the X,’s can be measured without prior 
knowledge of the mechanism, the evaluation of J* by difference 
hinges on the possibility of obtaining correct values of L,,. In 
general it is only then possible to evaluate the latter, and hence 
the former, when neither J* nor any of the L,, are functions of the 
X,’s. This is not the type of a priori assumption to be indiscrimi- 
nately applied to a system which may be expected to possess as 
large a number of feedback mechanisms as a living cell. Alterna- 
tively, if it may be assumed that L;,=9 for j 47, one can in one 
special case obtain a minimum value J*_,, from the consideration 
that L.. 2 0; if the directions are defined so that x, > 0, the smallest 
possible value of L,,X,is0. If the observed value J; < 0, Ji ee Jj. 
In widespread practice assumptions (ii) and (iii) are invoked in 
addition, allowing estimates when J.=0. None of these estimates 
are valid if L., 40, inasmuch as the restriction L;; 20 does not 
necessarily apply for L, ;j Zi. Since the possibility of ascertain- 
ing the occurrence of active” transport is dependent on an un- 
equivocal limiting value of the active flux term one is led to the 
conclusion that in the general case it is in principle impossible to 
ascertain, solely on the basis of flux and thermodynamic potential 
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measurements, whether a given species has been transported by a 
specific metabolic mechanism. Any statement as to the nature of 
the transport processes involved depends in its validity entirely 
on additional knowledge of the origin of the responsible thermo- 
dynamic forces. 


The author is indebted to Drs. F. M. Snell and C. S. Patlak for 
several valuable suggestions. 
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